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ABSTRACT 


The statistics of acoustic signals propagated to long ranges in the ocean 
are investigated in detail in this thesis. The phase random model of multipath 
propagation is extended to include finite bandwidth and/or modulated sources as 
well as multiple source configurations. The theoretical analyses include the | 
derivation of many new probability density functions for these new cases as well} 
as for the single narrowband source. | 


The probability density function for A, the time k of change of the levd 
in decibels for the single narrowband source is derived. : (A) depends only иро! 
у2, the single path mean Square phase rate, which can be A to certain oceal 
dynamical processes. The analysis of finite bandwidth and/or тот sources | 
reveals that the amplitude and amplitude rate densities (including P:(A)) are ` 

‚ independent of the finite bandwidth and modulation effects, but the M ON | 
the time rate-of- change of the multipath phase dis sepsitive to these effects. 
Thus, fitting P-(A) to histograms from data to find v^ is the preferred method 
for determining this important parameter. Bandwidth effects іп фсап be neglected 
however, if B<<2y where B is the signal bandwidth. The analysis atso reveals a 
potentially powerful technique for determining parameters of the modulation or 
bandwidth of a source from the received multipath signal. 


The analysis of multiple sources, applicable to noise problems, includes 
important approximations to densities which are intractable analytically, and 
would involve significant computer time to solve exactly. In addition to studies 
of the amplitude densities, significant progress has been made in solving for 
the amplitude rate densities and the joint densities of amplitude аы ampl itude 
rate. 


In addition to providing valuable confirmation of much of the theoretical 
analysis, a computer simulation of phase random multipath propagation also con- 
firms that for N > 4 paths phase random multipath conditions begin to closely 
approach the asymptotic conditions for N > œ. 


Data at 220Hz and 406Hz received by drifting sonobuoys in the Atlantic at 
approximately 300 km in range were analyzed. Values of y? obtained support an 
internal wave model for the relevant dynamical process. The modulation theory 
uncovered a heretofore unrecognized modulation in the data due to an error of 
the Doppler tracking system. Predictions of crossing rates including this modu- 
lation effect are in good agreement with the data. 


Other data at 15Hz and 33Hz propagated to ranges between 250 km to 450 km 
in the Pacific in which deliberate modulation was introduced, once again provide 
excellent confirmation of the theory. Measured values of v^ vary significantly 
from run to run and are not consistent with an internal wave model, indicating 
some other mechanism (i.e., tidal, rough scattering) must account for the fully 
saturated phase random nature of the data. The technique for determining modula- 
tion parameters was used, and for the 76 modulated runs analyzed, the average 
error in determining the actual bandwidth of the modulation of the source from 
the received multipath signal was 8%. 
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ABSTRACT 


The statistics of acoustic Signals propagated to long 
ranges in the ocean are investigated in detail in this 
thesis. The phase random model of multipath propagation 
1s extended to include finite bandwidth and/or modulated 
sources as well as multiple source configurations. The 
theoretical analyses include the derivation of many new 
probability density functions for these new cases as 
well as for the single narrowband source. 


The probability density function for A, the time 
D a Eor eange Of the level in decibels for the single 
narrowband source is derived. Р; (л) depends only upon e 
the single path mean square phase rate, which can be 
related to certain ocean dynamical processes. The analysis 
of finite bandwidth and/or modulated sources reveals that 
pnENnpistude and-ampiitude rate densities (including 
P'(A)) are independent of the finite bandwidth and 
modulation effects, but the density of the time 
eee Of Change of the multipath ря ф 15 sensitive to 
these effects. Thus, fitting P;(À) to histograms from 
ШИВ to find vs is the preferre а for determining 
this important parameter. Bandwidth effects in ó can be 
neglected, however, if B<< 2v where B is the signal 
bandwidth. The analysis also reveals a potentially 
powerful technique for determining parameters of the 





modulation or bandwidth of a source from the received 
multipath signal. 


The analysis of multiple sources, applicable to noise 
problems, includes important approximations to densities 
which are intractable analytically, and would involve 
significant computer time to solve exactly. In addition 
to studies of the amplitude densities, significant progress 
has been made in solving for the amplitude rate densities 
and the joint densities of amplitude and amplitude rate. 


In addition to providing valuable confirmation of 
much of the theoretical analysis, a computer simulation 
of phase random multipath propagation also confirms that 
for N > 4 paths phase random multipath conditions begin to 
closely approach the asymptotic conditions for N > >, 


Data at 220Hz and 406Hz received by drifting 
sonobuoys in the Atlantic at approximately 300 km in range 
were analyzed. Values of v* obtained support an internal 
wave model for the relevant dynamical process. The 
modulation theory uncovered a heretofore unrecognized 
modulation in the data due to an error of the Doppler 
tracking system. Predictions of crossing rates including 
this modulation effect are in good agreement with the data. 


Other data at 15Hz and 33Hz propagated to ranges 
between 250 km to 450 km in the Pacific in which deliberate 
modulation was introduced, once again provide excellent 
confirmation of the theory. Measured values of v* vary 
een icantly from run to run and are not consistent with 
an internal wave model, indicating some other mechanism 
E s5idal, rough scattering) must account for the fully 
Saturated phase random nature of the data. The technique 
for determining modulation parameters was used, and for 
the 76 modulated runs analyzed, the average error in 
determining the actual bandwidth of the modulation of the 
source from the received multipath signal was 8$. 
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coefficient of skew 
coefficient of excess 
ношен s constant (= .5772...) 
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parabolic cylinder function 


defined as the expected value of random 
variable x 


the exponential integral 


complete elliptic integral of the second 
kind 


defined as 10 1081) е = 4.34... 
Gauss's hypergeometric series 
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phase crossing rate 

single path phase for the nth path 
Bessel's function of order v 

the modified Bessel function of order v 
complete elliptic integral of the first kind 
number of sources 

number of sources with different intensity 
and different mean square phase rates at 
the receiver 

number of sources with equal intensity and 
equal mean square phase rates at the 


receiver 


level in dB's of the short time average 
mean-square pressure (A = 10 10g, 9x) 


time rate of change of A 
temporal modulation function 
time rate of change of M(t) 


characteristic function of random 
variable x 


two-dimensional characteristic function of 
the joint pdf of random variables X1 
and X2 


the confluent hypergeometric function 
[sometimes written $6(a,6,2)] 


number of groups of sources, each group with 
LE, ПАШЕ Cee) Sources per group 


number of independent propagation paths 
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coge 


the long time average mean square pressure 
at the receiver for a single source 
cu = 201 zum o) 


the single path mean-square phase rate 
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orean lity density function (pdf) of 
random variable x 


the joint pdf of random variables X1 
and X2 


the pdf of random variable y for N groups 


the Laplace transform of the pdf of 
random variable x 


Single path amplitude at the receiver 
short time average root-mean-square pressure 
time rate of change of p 


variance of phase in Gaussian phase 
modulation 


one half of the long time average 
mean-square pressure at the receiver 
(ci = Му^/2) 

Whittaker's function 


quadrature components of the signal envelope 


time rate of change of the quadrature 
components 


degenerate hypergeometric series in two 
variables 


short time average mean-square pressure 
Gane) 


time rate of change of x 


multipath phase 





y (x) 
w* (x) 


multipath phase rate 
Euler psi function [w(x) 


digamma function [Y' (x) 
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INTRODUCTION 


Acoustic signals propagated to long ranges in the 
Ocean, tens to hundreds to thousands of kilometers, via 
all modes of propagation including surface ducts, the 
deep sound channel, or sea-surface and ocean-bottom 
reflections exhibit fluctuations in amplitude and phase 
which are now recognized to be dominated by the multipath 
interference of the acoustic field. The statistics of 
these fluctuations as well as their relationships to the 
dynamics of the ocean has been one focus of recent research 
in understanding this important physical process. The 
optimum design of sonars (e.g., the receiver operating 
characteristics), underwater communications devices, and 
in fact any system which operates via acoustic 
transmission in the sea depends upon the knowledge of the 
Statistical behavior of these transmissions. 

The recognition of the dominance of the multipath 
Structure on the statistics, or the assertion that long 
range multipath acoustic propagation in the ocean can be 
modelled as a phase random process has been established 
only within the last ten years, although the phase random 
process or random walk problem has been under study since 
Rayleigh (1880) [1], and is one of the classical problems 


of mathematics and physics. Bergmann (1946) [2] was among 
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the first to speculate that observed fluctuations in 
signal intensity might be attributed to the interference 
of many paths summing in random phase. Dyer (1970) [3] 
formally applied the theory of a phase random process or 
random walk problem to long range acoustic multipath 
propagation in the ocean, and was the first to investigate 
the statistics of log transformed variables. Dyer also 
showed that even in the presence of scattering randomness 
multipath interference would dominate the statistics. 
This research in fact indicated a basic shift from the 
scattering models of earlier research which are more 
appropriate for high frequencies and short ranges when 
multipath effects are m. important. Dyer also proposed 
a model of distant shipping noise based upon the precepts 
of phase random multipath propagation and continued this 
research in a later paper (1973) [4], and most recently in 
Mikhalevsky and Dyer (1978) [5], results of the latter 
being included as part of this thesis. This model, 
appropriate to distant shipping noise, assumes the noise 
in a band is dominated by narrowband lines discrete in 
frequency. 

Mark (1972) [6] investigated the statistics of the 
multipath propagation of finite bandwidth signals. 


Employing a systems approach, he derived general expressions 
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for the mean and variance of the received signal energy 
in terms of the correlation function of the time varying 
impulse response of the medium. He also showed that the 
Gemma Or Erlang probability density function is often a 
good approximation to the real pdf of the received 
energy. Much earlier, Nakagami [7] had noted the 
utility of an appropriately transformed Erlang pdf in 
approximating the densities of received HF electromagnetic 
radiation undergoing rapid fading, indicating the 
broadness of scope of the phase random model and its 
general applicability. All these efforts, however, 
concentrated on the amplitude (or related quantities) of 
the signal and did not address the amplitude rate or 
phase rate of the signal. 

fonguet-Higgins (1975) [8] in connection with research 
on random sea surface waves (another phase random process) 
introduced the joint pdf's of amplitude, amplitude rate, 
phase, and phase rate as well as the marginal densities to 
the growing body of knowledge of phase random processes. 
It remained for Hamblen (1977) [9] to £ormally extend the 
phase random analysis to the multipath acoustic propagation 
process, incorporating the results of Longuet-Higgins [8] 
and also S.O. Rice [10] whose extensive research on noise 


statistics were also applicable to the long range acoustic 
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propagation problem.  Hamblen [9] also established the 
dependence of the probability density functions for a 
single narrowband acoustic source on the two parameters 
91, one-half the mean square pressure at the receiver, 
ЕИН 2/2. the single path mean are rate. He also 
verified the basic results with data from an ocean 
acoustic propagation experiment. 

Concurrent with the development of the phase random 
model of multipath acoustic propagation, much research 
was and is being conducted on another important aspect 
of the problem, namely to discover what ocean dynamic 
processes are the driving mechanisms and how parameters 
of these ocean dynamic processes are related to the 
parameters of 'the acoustic field, c1 and ví. Most 
notable perhaps is the recent research of Dyson, Munx, 
and Zetler (1976) [11] who have proposed a theoretical 
model relating the dynamics of internal waves in the ocean 
to the fluctuation of the acoustic field. However, this 
model appears to have serious limitations at low 
frequencies. In fact, little research has been reported 
on the low frequency cases. 

In the following paragraphs I will introduce the 
research reported on in this thesis. There are three 


basic areas in which significant progress has been made in 
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Mmmeerstanding the statistical behavior of multipath 
acoustic propagation: (1) the single narrowband source, 
(2) finite bandwidth and/or modulated sources, and 

(3) multiple sources or receptions. An extensive 

analysis of data from acoustic experiments in the ocean 

as well as a computer simulation of phase random multipath 
acoustic propagation not only increase our confidence in 
this new understanding, but reveal new information for low 
frequency signals about the driving mechanism of v^, the 


single path mean square phase rate. 


The Single Narrowband Source 


I derive for a single narrowband source the pdf's 
for the time rate of change of the sta (short time average) 
mean square pressure, ў, and the time rate of change ог 


the level in decibels, А, results which are unique to this 


^ 
- 


I 


metion only of v^. This result is of particular 


thesis.. The pdf for À is independent of of and is a 
importance as it affords a method of measuring v^ from 
ocean acoustic data without error due to uncertainties in 
the signal carrier. Included in this analysis are the 
joint densities Ex y and PA 12229 and their 


characteristic functions, also unique to this thesis. 


This analysis in fact completes the family of first order 
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and joint pdf's for a single narrowband source. 

E EHudix A contains a list of all the first order pdf's, 
their characteristic functions, means, and variances as 
found by me, and earlier as found by others, so that a 


complete set can be referred to. 


Finite Bandwidth and/or Modulated Sources 

Many acoustic signals of interest received in the 
ocean have bandwidths which are not narrow, and carriers 
that may not be stable or may be deliberately modulated. 
In applying the phase random model of multipath acoustic 
propagation, it is necessary to assume that the spectrum 
of the received signal is narrow, and that in homodyning 
the signal the spectral mean is zero [8,9]. Clearly, 
the signals mentioned in the beginning of this paragraph 
would violate these assumptions. I show, however, that 
the amplitude and amplitude rate statistics (including 
Р; (Л) are independent of finite bandwidth and/or carrier 
modulation effects. The multipath phase rate, 5, 15 
sensitive to these effects, and it is the pdf for Ф which 
must be modified. 

не; рег for 5 is in fact a function of v^, as well as 


the bandwidth, and/or parameters of the modulation. 


Bandwidth effects can be neglected when B<< 2v, which is 
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M@m@ererore a criterion for what is meant by "narrowband" 
in the ocean. Furthermore, this analysis reveals a 
method of separating and understanding the source induced 
modulation independent of the ocean induced modulation 
се versa. In addition to solving for the pdf's of 

Ф in the presence of bandwidth and/or modulation, I solve 
for the crossing rate statistics of phase for these cases 


as well. 


Multiple Sources or Receptions 

иешеша = 05 of the received signal amplitude and 
amplitude rate when there is multiple source structure 
depends on the exact nature of the received multi-source 
signal and the analysis performed by the receiver. І 
consider two basic cases. First, I assume the 
receptions (one per source) are disjoint in frequency in 
the analysis band and can be separated and summed 
incoherently (that is, without concern for phase). I 
assume, therefore, that each of the receptions/sources 
are independent, thus the analysis band should not include 
harmonics of a signal already in the band. This type 
of analysis is motivated by the structure of distant 
shipping noise [4,5]. For the second case, I consider 


the receptions to be at or so close in frequency that they 
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must be summed coherently. Note that the dividing line 
between case one and case two depends upon the resolution 
of the receiver. Most all of the pdf's I derive for the 
amplitude and particularly the amplitude rate for these 
cases are unique to this thesis. Where I have been 
unable to derive the exact pdf analytically or in which 
the exact solution is extremely time consuming to obtain 
even with the aid of a computer, I have in most cases 
found an approximation based upon Edgeworth's series. 
This approximation is shown to be excellent in the main 
lobe of the density but performance is degraded in the 
tails. 

I use the analysis of coherent sources to model the 
effect of ocean ambient noise on the pdf's for a single 
narrowband source. The pdf's are expressed in terms of 
the SNR (signal to noise ratio). In light of this 
analysis, I am able, as well, to extend with only slight 
modification all the frequency disjoint multiple source 
Solutions to include the cases when both coherent and 
joint source structure is present in the analysis band, 
or for noise which is continuously distributed over the 
passband. 

The reader is forewarned that the sections of this 


thesis on multiple sources (see Table of Contents) are 
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lengthy and in many cases quite tedious due to the 
complexity of the analysis. Table III is provided (page 143) 
to aid in understanding the organization of this material 
and to aid the reader in finding that analysis which 

is most pertinent to his problem. The analysis of the 
statistics of multiple sources which is not immediately 
motivated by any current ocean acoustics problem of 
interest is presented for completeness with potential for 
application ure problems, even perhaps in areas 
КООШ Се to acoustics. For example, the solutions for 
Eo , uSing the terminology of Queuing Theory, are in 
fact the pdf's for the interarrival times of cascaded 
Poisson processes which are unique to this thesis and, 

to my knowledge, not to be found in Queuing Theory 


literature. 


Computer Simulation and Data Analysis 


A computer simulation of phase random multipath 
propagation is developed to assist in and to check the 
theoretical analysis. The simulation demonstrates the 
independence of the exact nature of the pdf for 8 
(the single path phase) to the statistics when there are 
at least four propagation paths from the source to the 


receiver. Of importance is the confirmation of the 
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Ely accepted criterion that N > 4 paths results in 
an almost fully saturated phase random process, and pdf's 
Ehe limiting case as № ~ œ suffice. The simulation 
allows me in Appendix C to compare the theoretical 
pdf's, many which I derived, for N « 3 to the computer 
generated histograms with excellent results. 

I analyze data from two ocean acoustic experiments. 
One was conducted near Eleuthera in which two CW signals 
at 220Hz and 406Hz were transmitted approximately 300 km 
northeast towards Bermuda and received by drifting 
sonobuoys [12]. These data support the theoretical 
pdf's derived in Chapter l. Furthermore, the modulation 
theory uncovered a heretofore unrecognized modulation 
in the phase rate data due to errors in the sonobuoy 
tracking system that has dramatic results on the 
Statistics for the phase crossing rates. These data are 
also consistent with a model of phase random multipath 
propagation resulting from the interaction of the acoustic 
field with internal waves [11]. 

I also analyze data taken in the Pacific in 1973 
known collectively as the CASE experiment [13] in which 
CW signals at 15Hz and 33Hz were propagated to ranges 
varying from 250 km to 450 km. These data were 


deliberatelv modulated and the predictions derived 
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theoretically in Chapter 1 are once again confirmed. The 
analysis I used was capable of determining the bandwidth 
of the modulation from the received multipath signal 

with an average error of 8$. Where previously there had 
been some problems in the consistency of the data with the 
phase random model, the modulation theory successfully 
lays these problems to rest. These data, however, are 

not consistent with the internal wave model and this 
analysis has brought sharply into focus a shortcoming 

in our current understanding of the driving mechanisms of 
v^, particularly for low frequency signals. Correlations 
I have derived appear to support, though tenuously without 
additional arch, either a rough scattering or tidal 
mechanism to account for the fully saturated nature o£ 

the CASE data. 

In Appendix D, the effects of amplitude parameter 
variation are discussed. This analysis is aimed at 
uncovering the effects of temporal variations in the 
total energy of the signal during the observation period 


on the statistics of the received multipath signal. 
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CHAPTER 1 


ANALYTICAL TREATMENT 


E 5ingle Narrowband Source 


For a narrowband signal the phase-random model of 
multipath acoustic propagation predicts that the sta root- 
mean-square pressure, p, is a Rayleigh distributed random 
variable [3,9], its rate, §, is Gaussian [9], the multipath 
Басе, Ф, is uniform, and its rate, ф, is distributed 
ООА пс to a density first given by M.S. Lonqguet-Higgins 
ПЕН The transformation x = p^ yieids the density for 
the sta mean-square pressure which is exponential [3], 
and following Dyer [3] the transformation A = 10 log "X 
gives the density for the level in decibels which is 
Log-Rayleigh. To complete the family of first order 
densities for a single narrowband source, I have derived 


E densities for X апа A. 


IN Derivation of O and Р, (Л) 
For phase random multipath propagation the joint 
density of thé sta rms pressure p and its rate 9 is given 


Bey [8,9], 
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and N is the number of propagation paths, r is the single 
path amplitude, and u is the long-time average mean-square 
pressure. 

ШОО o we have у = 200. To find the pdf 
ferobability density function) for X I first solve for 
the cumulative distribution function (14 of x by 


integrating over the joint density of p and 6, 


ЕС Ша топ (1.1). 
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Eos 
9X 
result with Equations (1.1) and (1.3) differentiating 


The pdf for v will be given by Р 200. Combining this 


under the integral sign and integrating once, I obtain 


the result 








X 








= exp |- ‚ о-=сў<о (1.4) 
A 404 ^v 20,7 
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The absolute value results from the analysis for ү < 0. 
Equation (1.4) is the Laplace pdf [4]. Тһе variance of x 


is given by 


The pdf for А 15 solved in the same manner. First I 
find the pdf for y where y = lnp* and thus у = 26/р, 


therefore, 


o 22 
2 
• д • е е 
Р. = Aus P • , а а Р > S 
947) = 35 | | о, (P16) dedé y>0 (1.5) 
OO 


As before, I differentiate under the integral sign and 


perform the remaining integration to obtain 


E = ISE 2 р -о<ү< (1.6) 


d ју + |”? 
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Making the final transformation А = ey, where 
e = 10 logige = 4,34... (this notation will be adhered to 


E Docremarnder of the thesis), I obtain the final result 
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This density is in fact of the same form as Ps (6) (see 
Equation A8 in Appendix A). The second moment of this 
density is infinite. Physical insight of this phenomenon 
can be obtained if one visualizes the random walk problem. 
The amplitude of the vector is o and its phase d. When 
the amplitude of 9 goes to zero as ina deep fade it is 
Pr, seen that the phase can undergo very rapid changes. 
Likewise because A = Soa 15 also clear that A can 
assume very large values when p is small. Thus both À 
and ф are governed by Tae same form of pdf and their 
Variances are infinite. It is also noteworthy that 


while the pdf's for f and X depend upon both c, and ve 


1 
Equation (1.7) depends only on v. 


ШОО 2 Joint Densities of Amplitude and Amplitude Rate 
ammerossing Rate Statistics 


In order to complete the family of joint densities 
of amplitude and amplitude rate for a single source, I 
have derived Р, y XX) and Р, 104,49. Thıs analysıs also 
provides an alternative method for deriving Pe (Xx) and 


BELA) to check the calculations of the previous section. 


In addition, I derive the two-dimensional characteristic 
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Benetzions 4] which will be of significant importance when 


I investigate multiple sources. 


For Х = p? we have р = XE ana p = 5 ў ELE 
order to solve for Р, (ХХ) I make a two dimensional 
, 
transformation on ЕЗ a (PÒ) (Equation 1.1) as follows: 
Iı 


nmel =172 


22 | 7 ХХ 
Р . х) = P . оаа й 
y AO | 6,5 (9,5) раб (1.8) 
О = оо 


By differentiating under the integrals in Equation (1.8) 


ENcb5btain 
i X x? 
E. Rx) e à expí- = 5 zu 
, 
401 уу 2ту 291 8x0, V 
х > 0, and -œ<%<% (1.9) 


O dimensional characteristic function or two 


dimensional Fourier transform can be defined as 
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куха. | | | хх) Хака) е dx,dx, (1.10) 
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Performing the double integration yields 


M. (0,0) ч -------------- (T2) 


Recalling that y = lny, we have y = eY and у = ўе?. 
Thus, 
m yet 
. 9 : : 
P o , = —— Р S , а а 10013 
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As before I differentiate under the integrals, and make the 


final transformations A = ey and Å = eý to obtain, 
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The two dimensional characteristic function of РА д ОЛА) 
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(Equation 1.14) is 
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where E, (x) is the modified Bessel function of order 2. 
Performing the integrals over x and A in Equations 
(1.9) and (1.14) respectively yields the marginal densities 
of x and Ä which are given by Equations (1.4) and: (1.7) 
respectively as expected. 
Following Rice [10] and using Equation (1.14) I 
derive the mean axis crossing rate for the amplitude of 
the signal expressed in decibels. 


BL. еі заан | 
G(A,) -2| ір, 10 -,Маі ғ ELA A= Ag! (1.16) 
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If we transform Equation (1.17) to determine the axis 


crossing rate for 22 = exp (4/28) 
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Equation (1.18) was previously obtained by Dyer and 


Shepard [15], and Hamblen [9]. 


1.2 Multiple Sources 
In this section the statistics of the amplitude and 


amplitude rate variables are analyzed when the source 
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structure consists of many independent acoustic generators 
distributed at various frequencies across the bandwidth 

of the receiver. It is assumed that propagation from each 
source is phase random and multipath. An initial assump- 
tion is also made that each source radiates at a different 
frequency. This problem is applicable to the case of deep 
ocean ambient noise due to distant shipping [4]. For 
this case the spectrum of ship radiation is assumed to be 
dominated by lines disjoint in frequency, so that each 
ship contributes as many independent sources as 

there are lines within the observational bandwidth. 

(Note, the observational bandwidth cannot be so large as 
to include harmonics.) Henceforth line and source will 

be used interchangeably. Furthermore, the analysis 
assumes that the Fourier components of the received 

signal over the observational bandwidth are squared and 
summed. Thus, by Parseval's theorem it is the square 


Of the signal amplitude, in this case у = p? which is 


L 
ES Xeotal Е nél Xn 


As noted in Reference [4], this model breaks down when 


summed for each source, Tori, sources. 
sources cannot be separated in frequency, and then the 
model must be modified to include the effect of two or 
more sources which may be radiating at the same frequency. 


Hits coherent problem is treated in Section 1.2.3. 
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1.2.1 Amplitude Densities 

Dyer [4] first investigated the amplitude statistics 
for the multicomponent case and developed the framework 
for the problem which will be followed here. Three 
categories are defined: (a) all line components, L in 
number, arrive with the same long-time average intensity; 
(b) all line components, M in number, arrive with different 
intensities; and (c) N groups arrive, each with L; equal 
intensity components. As an example, Case (c) was 
applied in Reference [4] to noise as might be measured 
at low frequencies in deep water near Bermuda. 

As pointed out in Reference [4], and as will be 
shown in Section 1.2.1.3, use of the models describing 
line component noise often entails considerable 
computational tedium. This complexity often motivates 
adoption of approximate methods, which will be discussed 
in Section 1.2.1.4. Cases (a) and (b) are treated in 
Reference [4]. Exact solutions for N=2 and N=3 of 
ece (c) are derived in Sections 1.2.1.1 and 1.2.1.2 
respectively. It is true that the exact solutions of 
N22 and N23 may be of little practical value as most 
‘cases of interest will probably contain many more than 
three groups. The analysis is performed, however, to 


ENG the theoretical framework of this problem, perhaps 
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enabling someone else to solve it exactly in all 
generality, and more importantly in relation to this 
thesis research these exact solutions will provide a check 
for the approximate solution which is presented in 


BEctron 1.2.1.4. 


ID UI Exact Solution for N=2, Case (c) 
Consider L line components, each with the same 


long-time average mean square pressure y. The 


probability of the sta mean square x is [3 4] 


au 1-1.-ах 
Р (= E , Х>0,а>0 ; 
XLE (L-1)! 
Ecco а | (151.9) 
where a = l/u, the mean Hy = Lue and the variance 
L 
т = ius = Uu /L. Equation (1.19) is the Erlang or 
L L 


gamma pdf (see Appendix B). 

Let X4 and X» be distributed according to 
Esnation (1.19) but with different ні and different L; 
= 1,2). Further, let X = X1 +X. Then 


X 
P y VO = f Ру (5) Р.(Х-5)95. If we make the change of 
2 О 
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variables t = xt, 
Б Lo = OX 1 
а а е L-+L.-1 Ше Бе | 
P(X) Киа сте н 2 | eL Te) 2 е®°Х(а2-а1) + 
2 T(14) l'(L5) 0 | 


ALIZO) 
The integral is the confluent hypergeometric function 
defined in its integral form [16,17]: 
1 
NON --—l _ | e2* pth yg STO} ae (1.21) 
Ic coo) 0 


ЖЕНЕ 15 a well tabulated function [18 ,19,20] alternatively 


defined by Kummers series (absolutely convergent) : 
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where (a). = а(а+1) (а+2)... (а+п-1), (а). = 


A се (1.21) to Equation (1.20) I obtain 
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the desired result: 


a, “La Vay beth 197 X22 
BENED — —— — — — M U45424,*25^;Xx(8578,)] (1.23) 
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Ef a, = a, then as expected, Equation (1.23) reduces to 
Eeuation (1.19) with L = Li + La y upon noting from 
ecuacion (1.22) that M(a,ô,0) = 1. 

To obtain the density of the level in dB I make the 


now familiar transformation A = = lny in Equation (1.23): 


а) Ла, 2 A 
P AU) = 1 [L, ,L,+L,, (a,-a,)e' FE), 
2 eT(L, +L,) 


exp > (L + L) - а. QOIS (2 0) 





The statistics for A are: 


н, = е Јаџ, - y + S41 (I4 * L5) E Q1 (44,85, L4 ;D5)] 
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| а Ly © (L.) a.-a Nu+n-1 k-1 
T 1 m 27751 1 
Q,(a] ray, L],L,) = 2) ) = Е у у x 
2 n=2 ` 2 к=Мм+] o=M 
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where M = Li ~ — ӘНЕ = Euler s constant = .5772... 


Note that if а, = а„, Equations (1.25) reduce to 
Equations (4) of Reference [4] as expected. 


For non-integer values of Li and Lo, Equations (1,25) 


can be expressed as follows: 








а n о (Li). a,-a, S 
u,/e = inu, + |— I [y (L¿+L, +n)] 





а, =0 ni a, 
Gb. 26) 
L n 
a 1 о CES anra 
2 2 EI i er 2 1 { _ 2 
(04 * u1)/& 2 Pru => | [P(L,+L,+n) 1naj] 








E y? (L,+L,+n) | 


where 
W(x) is the Euler Psi function [21] 
and where 


Wood = 22 #00. 





ge 


Equation (1.23) can also be expressed in terms of the 
incomplete gamma function [22,23] which is a special case 


of the confluent hypergeometric function. 


ШОО 1.2 Exact Solution for N23, Case (c) 

I now consider three groups with L; er = 1,2, s)E sources 
in each group. Let 4X * X4 * X4 * X4 = 9X * X4 where the 
Xi 2,3) axe distributed according to Equation (1.19) 
(with different means and L,) and „x is distributed 


Beoardına to Equation (1.23). Then: 


X 
P - Р Р = уа 
y VO | К) A сас 


3 о 2 


As before, let z = xt and 


Li L5 L, L,+L,+L,-1 -а-х 
a) a, a, X e 
P x =  — _Á s I (2070) 
3 P (Li tL) T (L3) 


and 





| 
Dee = ty(a,-a,) 
Е B 3 ger 2 E Edu. 
Т = | (l-t) Ё е M [Lj ,LitL5, 
0 


tx (a,-a 


ees 
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(1228) 


Equation (1.28) is evaluated by expanding the confluent 


hypergeometric junction in its series form (Equation 1.22), 


and integrating term by term. When I combine this 


result with Equation (1.27) I obtain the result 


L G PL LE S -a4X 


Е l- 2 а 3 X т ду = a = (а E yn 
NE 2 3 PN: 
EM ————— —— MELLE 
3 P (Ly +L,+L3) n=0 
o M[L, +L,+n, L} +L,+L +n, x(a,-a,)] 
with 
E lE и 
Жаз 2 2 2 
x A MU LE 


(L 
+L 


en 
(Т "m 


E n 


S 


Deriving numerical results from Equation (1.29) is 


Straightforward by rewriting the summation 





nzo PP 
where 
| x(a,-a,) L,+n-1 
Sn De “а-і) “o”? Zn 
Irre 
and [18], 
b(lrb-2z) psc 
DELE m _]ı + —— m 6: (15582) 
az az 


in these I define m, = M (a+n, b+n, z) and, a = +L 


L 
IN 2.2 
b = L,+L,+L,, апа 2 - x (a4 7 a4). 
From Equation (1.29) the density of the level in dB 


is obtained as before: 


Р (А) = = eYp  (eY) (1.33) 
3 3% у = Me 


The statistics of the level are: 
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= 
= 
І 


= € [13nu4 7 Y * $4 (L4 *L5*L,) = R, (à4,a5,a4,L4, L5, L4)] 


2 
E = € Е S,(L]+L,+L,) R1 (a,,a5,a4, L4 ; L4; b4) 








* КО 12,15) | (1.34) 
where: 
L L D 
1 2 
a a o  (L,) a.-a 
1 2 SO 2 tc E 
R. (a,,a,,a,,L,,L,,L,) = |— т ) zz 
a 25 1 23 = = p=0 p! = 
© (L,+L.+p) Ja) penel 
dh о 1 
BC n! а; kz0 Lj*L4*b4*X 


ког реп > 0, and 
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a a со L алта 
1 2 1 теке ДІ 
Bar ,a,.,2,,L.,L,.,L,) = 2 — = ў E 
НЕ 2" 3 1 23 = = о=0 р! | a, 
со = =] = 
05 (L,+L)+p), |a4-a, ud 1 Y 1 
n=0 п: 55 ігі mag UT 
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Ern» Ii, in which # = L] +L, +L}. 
Note that for both N= 2 and N=3 one can judiciously 


choose which а; ,1; will be designated ayı li a,b or 
azr D, so as to insure the fastest convergence of the terms 


involving infinite series. 


uu 3 SOlution for Arbitrary N, Case (c) 

In this section I will derive an expression for the 
pdf of the sta mean square pressure у, for an arbitrary 
number of groups N, using Laplace transform techniques. 


дес 


М 
= 2 Ka (1.35) 


where NX denotes the random variable whose pdf we seek 
and the X, are Erlang distributed random variables, 
tributed according to Equation (1.19) with 


arbitrary order La’ The Laplace transform of у 15 


L 


n 
a 
А n 
ex i : + а ‚ые ксл (1.36) 


where m 1/u,- From the properties of Laplace transforms 








Zoe 


and independent random variables 


О N а, La 
(5) = I (————) 2555 70 
NA n=1 ? i an 
End thus, 
je 
N a L 
1. n n Sx 
Р (x) = = | П ——— | = а5 (2258) 
NX 27) іе на 


BEhosngh Equation (1.38) is attractively compact, it is 
computationally бе лз for most cases of interest, i.e., 
large N and large L as you might expect with many ships 
and many lines. To carry the analysis a bit further, I 


make a partial fraction expansion of Equation (1.37), 


L 
y C hk 


N а N 


Nh. п=1 К 


It is clear from this expression that P x o9 is a linear 
N 
sum of weighted Erlang distributions, and the problem is 





EE 


reduced to finding the coefficients Car” I take the 
inverse transform of Equation (1.39) by applying the 


calculus of residues to obtain 


k-1 


за k 
Ck * ent „ar ta) 9.0 (1.40) 
s N "mA 
m 
and finally, 
СО = T "n y + Cu e (1.41) 
x) = а = | : 
NX nel ® деј кеј с 


The generally intractable nature of Equation (1.40) leads 
to the introduction of approximate methods discussed in 


the following section. 


1.2.1.4 Edgeworth s Series Approximation, Case (с) 


As before, 


where the X, are independent random variables distributed 





 бҮсіпс to Equation (1.19) and in general with 

different long term average means Hi and different 

number of line components 1.. Although ka cannot be 
solved exactly for N > 3 without considerable computational 
tedium, it can be approximated in the main lobe, and as 


we shall see, quite accurately and easily by an 


Edgeworth's series [М]: 


Em 1 l (3) M. (4) 
BE DO Hic joss) = = 7-2 (5) + 77 ¥2 (5) 
N NX 
610 2206) (ку - 1 25 105) (2) 
6! S 51 g? 
X 
- 27 vvzU (Е) 
= 455 iz (5) + ... terms in higher order 
moments i-o. 


where 








SI 


й E 
ШО уу Be gs 
2000 cm А е , 
С x (2т) 172 
М 
ED (cy - n derivative of the Gaussian pdf, 
2m K, 
ү. = coefficient of skew = —A = — (1.43) 
Ge Ge 
X X 
"y Ка 
ШЕ = coefficient of excess = —-3= — (5.4412 
© ap" 
X X 
"uk z nen central moment, and 
Ka = nen order cumulant or semi-invariant. 


The point-by-point error in the approximation of 
Equation (1.42) is of the same order as the first term 
neglected [4]. 

A computationally more efficient expression of 
Equation (1.42) can be derived by actually taking the 
derivatives of Z(£) as indicated with respect to £ and 
collect terms in powers of §. I have used only the first 
four terms in Equation (1.42) and as will he demonstrated 


later, this will be entirely adequate for practical 








aie 


computations. 
зү Y 
E e. SE 
N NX 


ПАО Бе (1.45) 


It should be noted that because Equation (1.45) is an 
approximation, negative values may be obtained for some 
regions in the tails. In fact, the Edgeworth's series 
performs best in the main lobe of the density and worse 
in the tails. The mean and variance is respectively 
(evoking the properties of sums of independent random 


variables f4)): 


N 
D y Lou. (1.46) 
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For the gamma ра: [14]: 

к, = wu LT (n) (1.48) 
Furthermore, for a sum of independent random variables [14]: 
K, + 3X, + +--+ + yK (1.49) 


MESS Equations (1.46) through (1.49) to obtain the 


coefficients of skew and excess for P . (x): 


N 
1 : 
Y = == ОШОН: , (1.50) 
= g? i=l a 
x 
1 | 
= o LM (су, 
E oy i=l = 


= and ү. аге zero for the Gaussian pdf and attain maximum 
values of 2 and 6 respectively for the exponential pact 
(which is the case of N=1 and L=]). 

Before examining the nature of ve and Ya more closely, 
let me return to Equation (1.42) and make the log 


transformation to obtain the Edgeworth's approximation for 





сое 


the level in dB: 











у y (801106 
pm) e e Ian - 2" ma 2) qq) 4 2 (п) 
E Р ! ! 
N 
Т2) 
or alternatively, 
зү. ЗҮ Y о ен У 
1 A/e EN S_ S SS 
EU ЕЕЕ ІІ ais Sa N 
X 
N 
Y MM Ye 10ү? 
- 3T n! (3T - Zn Es D ea (1.53) 
where 
nc Босо e) = u 
NA N^ 


As a result of the fact that the transformation 
А = 10 log, ox is nonlinear, the statistics of A cannot be 


found without exact knowledge of P 69: 
М 








ds 


(1,52) 


Using Equation (1.41) in the previous section and inter- 


changing the order of summation and integration I obtain: 


N L у Ж Chk j K-l _-anX 
ИЕ Сс ПШ, ат) ett | mex a 0 
NA п=1 n ЕЗІ 21 (k-1)! | 
(2252 
L 
N N n С 
I k-l _-a 
EN = (ce? п an) ) ) zu ems an) e nXay 
м^ м^ n=1 ™ n=l k=l (k=l): 
Evaluating Equations (1.55) yields [24]: 
N L г Cok 
ШІ - с Па"); —— [v(k) - 1n a,] 
м^ п=1 ™ n=l k=l a“ д 
(156) 
L 
N N noe 
ДЕ“, = de T arn) } | (too c + (x) | 
N N п=] п=] k=la, ) 
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Mmaomgh Equations (1.56) are exact, they are predicated on 
Ehe knowledge of the Cag S Avene Equation (1.40), 
which in turn motivates an approximate solution. 
Unfortunately, direct integration of the Edgeworth's 
series is not possible because term by term the integrals 
diverge. An approximation can be made, however, using 
MN? st term only, that is, to integrate the log 


transformed Gaussian. Thus, 





1257) 


S (х-н) 
NU Qo | amate - Ts 
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The integrals in Equations (1.57) are not straightforward 
but can be evaluated as follows. First, the square in the 
exponential is expanded, and the constant term removed 

С under the integral sign. Second, the change of 


variables x - nae is made: 





Zen 


u 
=372 со EE Т 
2 vy ce X | n inyu exp(- — 4 —& /u)du 
N 2c у2т Dias s 
X О Х Х 
(1558) 
12 
2. a u 
tu, „52 x | u “ОЛ к хе dH 
N N 2c 2T да ое 
Х О x Х 
Third, the radical is removed from the ln and the 
exponential of Yu is expanded in its power series. 
Fourth, the order of integration and summation are 
interchanged and finally I have, 
12 
^, ce 29X ENTE n қ (2т-2) u 
"A 2 —-- = E | u Inu exp (-——)du 
N 4с У/2т m=0 5 0) 5 
Х О x 
(1759) 
12 
- 2. 


NEU 20 > u, um (=m-=) 
N 20, 


The integral in Equations (1.59) can be evaluated [24] and 


the final result is obtained: 
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202 ә 1 У2ы 4m 

E = ге X ES = гал» Әуе + 1 207} 

N Дуп ш=0 ^" 

апа (1.60a) 

12 
2 ос 99 You m 
ШО - 5 Ку = rn 1) 
N N 8T m=0 ^" X 


: TEIE 2 1n207]^ = y Gard } 


Although the Edgeworth's series will provide an accurate 
approximation to the density we seek, the best approxima- 
tion for the moments of A are those of the log transformed 
Gaussian given by the statistics of Equations (1.60). 
Though not immediately derivable from Equations (1.60) 
when o is small, Dyer [4] obtains, by taking an expansion 
around the peak of the density, 
HA = (Inu, - (0/24 le 

ІЗ 20110 


2 x 2 2 2 
O = (0 
К ( "D € 





Returning to the Edgeworth's series, and to gain more 
insight into the way the coefficients of skew and excess 
EY. and pe behave, depending upon the number of groups N, 
and number of line components in each group Li, I have 
plotted Ys VS L пса l for N Tand N> 3. In like 
manner, Ya is plotted in Figure 2. Each was constructed 
as follows: Each group N has the same number of line 
components L, and the long term average mean у; of each 


Bzroup is quantitized in 3dB steps, i.e.: 


1 1-1 
ша = B relative units 


For the case N= 1, Equations (1.50) and (1.51) yield: 


E -1/2 
= 2L 

— -1 
Еке 6L 


En the limit as N-, Equations (1.50) and (1.51) 


 Пуетде to: 


= [1.48]1, 2/2 
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Fig.2 The coefficient of excess, yp, for one group (N=!) 
and three or more groups (N23) as a function 
of L , the number of line components in each group. 
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The convergence is very rapid. As evidence, the exact 
value for N=3 is only about 10% greater than the limiting 
value. 

Prom Figures l and 2 it can be seen that when N= 1, 
the coefficient of skew and excess closely approximate 
ШОО ап values when L> 6. When N>3 they closely 
E soximate Gaussian values for L» 3. More justification 
for this conclusion is given in the next several 
B uns. 

БЕ іптспег illustration of significance of the value 
of D. and A I consider the case N=1 and L=2. For 


по теслеу, I let u=1. Then, 


ex 


Р, (Х) = уе (1.61) 
= 2 
E. 
g? = 2 
X 
2 /2 
> 


Equation (1.61) is plotted with the Edgeworth's series, 
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Equation (1.45), and its corresponding Gaussian in 
Figure 3. Note that the large values of У апа ШЕ would 
predict the Gaussian to be a poor approximation, and so 
it is. The Edgeworth's series, on the other hand, is 
quite good except at the tails. 

mung the log transformation of Equation (1.61) for 
the same simple case (N=1 and L=2), I have the density 


of the level in dB: 


P (4) = = ехр{21 = exp (ae) } (1.62) 


Equation (1.62) is plotted with the Edgeworth's series as 
given by Equation (1.53) in Figure 4. Again, the values 
Eor Di. and Y, Suggest the log transform of the Gaussian 

to be a poor approximation, and indeed it is. But the 
Edgeworth's series for the density of the level is 
remarkably close. Thus, I conclude than when 2 апа Ye 
are large, the Gaussian is not a useful approximation, 

but the Edgeworth is, especially when dealing with the pdf 
of the level. Also, it is interesting to compare 
estimates of the mean level. Equation (1.603) gives 


Hp = 2.28dB while the exact value given by Equation (1.25), 
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Fig.3 Probability density for the case N=!, L=2,and =l. 
Ys -/2 апа Ye? 3. The exact density is shown 
with its Edgeworth approximation and its 


corresponding Gaussian. 
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Fig.4 Probability density of the level for the case plotted in Fig. 3. 
Shown with the exact density are its transformed 
Edgeworth approximation and its transformed Gaussian. 
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Setting Ly = L, = ау = а. 


As a second example, I consider Case (c) of Figure 2 


= 1, 15 и = пао 


of Reference [4], with the specification N = 2, a, = DO 


a, = 956, and Ly = L, = 2. For this we have (а= 1.14 
апа E = 2.04. The log transformed Edgeworth approximation, 
Equation (1.53), is plotted with the exact density in 
Figure 5. Here, Equation (1.60a) gives pues 9.16dB while 
the exact value given by Equation (1.25) is Ly = goo B 
The log transform of the Gaussian is not a good approxima- 
tion, but here again the Edgeworth and the exact density 
are for all practical purposes identical. 

For a final example I examine a three group problem 


(N= 3) using the results of Section 1.2.1.2 to compare 


with the Edgeworth's approximation and the Gaussian. 


BN this example, ај = i a, = Ža а. = 87 Là = 2, L, = 4, 
апа L4 ae Using Equatsons (1.49) and (1.50), I find 
Е = .86 and U IN MEC uUation (1.53). and the exact 


density of the level as given by Equation (1.33) are 
plotted in Figure 6 with the transformed Gaussian. Once 
again, Equation (1.53) is exact for all practical purposes 
while the Gaussian assumption will result in some error. 
Thus, from our examples for Ys and А МОБ too large, 
EL 0s 1.5 and 3 respectively, the first four terms of 


the log transformed Edgeworth series adecuately represents 
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Fig.9 Probability density of the level for the case N=2, 
LL: | p, * 0/5, ond E57 9/3, Ys: 1.14 and 
VIE 2.04 . Shown with the exact density are its 
transformed Edgeworth approximation and its 
transformed Gaussian. 
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Fig. 6 Probability density of the level for the case N=3, 
Piece Dr L,:16 zu eu ie опа 
3 1/8. y¿=.86 and ye= 1.28. Shown with the 
exact density are its transformed Edgeworth 
approximation and its transformed Gaussian. 
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the levels, while the log transformed Gaussian requires 
Т. 8; A % 0 for it to be adequate. 

I now apply the foregoing analysis to the N= 15 
problem of Reference [4]. This problem involves 
estimates of ocean traffic in the пећке clans 
Components from each of the ships are grouped with 3dB 
quantitization in N= l5 steps, with components within each 
group numbering as high as L; = 57, as is shown in the 
Ш ЕЗЕ three columns of Table I. With the use of 
macions (1.50) and (1.51) I find us отапа 


у_ = 1.743. In Table I, the contributions of each of the 


е 


groups to the mean, variance, and coefficients of skew 


rand excess are also tabulated. As is clear from the 


table, groups 9-15 contribute very little to the overall 
ty, since the variance, skew, and excess do not 

change (to within three decimal places) beyond N= 8. 
Reduction to an eight group problem, however, is not 

much of an improvement over the 15 group one. But, 
Equations (1.46), (1.47), (1.50), and (1.51) are very 
simple for any N, and in comparing ye and ШЕ for this 

case with those of our previous examples, we can expect 

the density of the mean square pressure to be closely 
MeproOximated by Equation (1.45). Similarly, Equation (1.53) 


should yield the density of the level in GB, which for all 
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Medetical purposes would be exact. Equation (1.45) for 
this case is plotted along with the Gaussian in Figure 7, 
and Equation (1.53) is plotted with the transformed 
Gaussian in Figure 8. The Gaussian pdf is seen to depart 
significantly from the Edgeworth pdf, even for a case 
having a very large number of line components such as 
may be appropriate to an actual oceanic situation. 
However, with reference to the examples above and 

Figures 3-6, we are justified in expecting that the 
Edgeworth pdf is virtually exact. 

| As noted earlier, an analytical expression has not 
yet been derived for the statistics of the log transform 
for N > 3. However, I can estimate the mean from 

EE. 8 to be Пе by a numerical integration of 


B uation 1.53) 


Un % 7.2ав 

Тһе Gaussian assumption used in Reference [4] led to a 
slightly higher value for the mean (% 7.5dB) as a glance 

at Figure 8 would explain. Equation (1.603 gives 

Ba % 7.51dB affirming the result obtained by Dyer [4]. The 
standard deviation as computed in Reference [4] of 1.2dB 
appears reasonable, again by inspection of Figure 8, while 


Fauation (1.60a) vields с, = 1.40dB. 
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Fig.” The Gaussian and Edgeworth approximations of the 
probability density of the short-time average mean 
square pressure of the noise that may be sensed 
in deep water near Bermuda in winter , for a |/3-oct 
band at 60 Hz and an omni -directional hydrophone. 
The ships have been grouped in 3dB steps in 
fifteen groups, with overall mean , variance , and 
coefficients of skew and excess as specified in 
Table I. 
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Fig.8 The transformed Gaussian and the transformed 
Edgeworth approximation of the probability density 
of the level for the case plotted in Fig. 7. 
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We may conclude that when the number of groups which 
add significantly to the overall density is greater than 
B, then use of the Edgeworth ЗИРЕ approximation сап be 
quite valuable in estimating the pdf of the received 
signal. While it may be tempting to use a Gaussian pdf 
for values of Ves l and Yo 2, Significant differences 
in the pdf for the mean square pressure or the level 
should dissuade us from this course. It is true that 
the mean and standard deviation of the levels are less 
sensitive to the differences between the Edgeworth and 
the Gaussian approximations, but the Edgeworth is not 
much more difficult to use and is thus to be recommended. 

The Edgeworth's series approximation therefore 
provides an easily implemented method of (1) estimating 
the statistics of the level in dB for even the most 
complicated realization of Case (c), and (2) revealing 
the extent to which the Gaussian assumption is a valid one. 
This is particularly valuable because one cannot merely 
assume on the basis that N is large that the Gaussian 
assumption will be a valid one. The critical factor is 
the amount of energy in each group. If one or two groups 
contain most of the energy, then the governing density 
will be significantly different from the Gaussian and, in 


fact, will more closely resembie the density that would be 
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associated with the most energetic group. 

If the tails of the pdf are the primary regions of 
interest, then exact computer solution may be a viable 
alternative for a complicated realization of Case (c). 
It should be noted, however, that more terms in the 
Edgeworth's series can be m to obtain any arbitrary 
accuracy desired, or alternative methods such as the 
Chernoff bound or "tilting" the density can be applied 


[25]. These methods will not be discussed in this work. 


1.2.2 Amplitude Rate Densities 


In this section of the thesis I will derive the 
pdf's of 6, ў, and A for the multiple source case. 
These rate variables are dependent variables with respect 
to the amplitude variables except for certain special 
cases. This fact introduces a great deal more 
complexity than has been encountered up to now. Solution 
for the joint densities of amplitude and amplitude rate 
are generally required before the marginal rate densities 
themselves can be found. Solution of the rate density 
for one variable does not lead by simple transformation 
to the solution for the rate densities for the other two 


variables as was the case for the amplitude. 


2 


In general, the rate varıables depend upon both 01; 
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| пеазџке ОБ Еће intensity, and v*, the single path mean 
square phase rate. This implies that there are a plethora 
of different possible combinations of multiple sources 
with different or same сі and/or different or same у7. 
The breakdown into the same cases employed in Section 1.2.1 
will be followed here, though their definitions must be 
expanded to include v*. In Section 1.2.2.1 I examine 
Case (a), the case of multiple sources or lines in which 
91 and v^ are the same for all source/receiver pairs. 
This would apply, for example, to noise which is flat 
across the passband of the receiver from a small 
geographical area or sector. In Section 1.2.2.2 I 
investigate various special cases when the 01'5 and the 

2 


v^ may be different for each source/receiver pair, with 


more general applicability, including Cases (b) and (c). 


ШРЕК Multiple Components ОЁ Equal eee tone 
and Equal осше Path Phase Rate 


will first solve for E (x) which is the pdf for 


LE 
X from L sources or components of equal oF and v^. From 
L L 
E ral “= E Nor it is clear that ee AT Dos: а. Thus 


mee this, and only this variable, we can employ all the 
very nice properties of sums of independent random 


Variables. For the pdf of XLE 





ge 


Yi Ya: X 
E. = | | | me en) No, 
m5 
Dr 
И 2 е 
- 1 р lox 
P. (у) = == | ————— << do (1764) 
XLE 21 (b* + "E 


- ОО 


by making use of M, (w) (Equation A6), where b = 1/201v. 


Por x > 0, therefore, 


2L L-1 iox 
: b JE d = 
ПИ 2 e (1.65) 
XLE 27 с aa әш 





w= ib 


Applying contour integration around the lower half plane 
Бог x< 0 yields, as expected, the identical result given 
by Equation (1.65) because x is symmetric. Applying 


Semetions (1.63) and (1.65) for L = 1-6, I obtain 


А b А 
Р. (X) = = exp(-b|x|) 
Х1Е e 





below) 


X6E 


(X) 


1 
(1-1) 12° кеј (k-1)!(L-k)!2 


zone 
ш r 2 | • 

mtb Ix| * 5) exp(-b| XD) 

{ь? | |2 + 3b? |X| - 3bJexp(-b| X |) 


(b*|x|? -* бр ||: + 15b? |X| + 15bJexp(-b| Xl) 


T (b^|x|*-10ob'Ix|?*-45b!lxX|? * 105p? |&l 


105b) ехр(-Ь|%|) 


| 


sr А ин Бе e aa e 


9455: || + 9455) ехр(-Ь|Х|) 


.(1.662-£) 


Upon inspection of Equations (1.66) the general form of 


the density for arbitrary L emerges (which will be proven 


5 L 
Au ХІ (L+k-2) ! 


mbar) 
xci |X! 


b 


20.07) 
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ETearly, u. = 0, and in general 
ÄLE 
en и 1 ч MEME .L-k+n _-byY .. 
Be: = L-1 k-1 | X E ах 
(1-1) !2 k=1 (k-1)! (L-k)!2 A 
GS) 
СП, for n odd is 0, and for n even: 
un 1 ern 
A O RAR 59) 
(1-1) 12 b k=l (k-1)!(L-k)!2 
Using the independence property (14 I get 
„2 к 50 E A 
E [Xr £f] = P SLO] V (1.70) 


This result is identical to Equation (1.69) when n=2. 


A much more elegant solution for 29 (X) will now be 
LE 
developed. The approach yields the complete solution for 
1 equal v? case including P. (5), D (A), and 
PLE LE 
all the joint densities of amplitude and amplitude rate. 


the equal ос 





Lo 


L 
ШОО п with Equation (1.12). Because Хүр ” ) Yo 
i ГЕ 
%-- ә) X, the x 's are independent of 
LE E] n n 
each other, and the XQ S are likewise independent of each 


other the joint density Е У (x,x) can be expressed as 
OTE 
follows: 





L 
| Zi VEN SION 
p ; (x,%) = 3T — | | ————á (ШЕЛІ) 
Xr KLE DE [ш + 1 (5+ 2020202) ] 
-— = E 1 
which is just the inverse transform of the = power of 


Eauation (1.12). I use the calculus of residues and 


perform the integration over w first, which has a pole of 











erder L: 
L со 
L-1 he 
с 2m: d =IWX zu, 
EN (Х.Х) = = ————— | = le | e g 
SALE - A e 
DES d M pc aD 
= 207 21У72 (1:72) 


and 





L со Ẹ 
L-1 - (—r + 
Е +: 2с 
Xr ALE mE а) 


Performing the final integration I obtain 


result, the joint density of the sta mean 


y and its rate Х for L sources with equal 
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29100) х - ісў 
ас 


O a) 


the desired 
square pressure 


91 and equal v: 


L-3 , 

X 
E ox - Sp oz 
SCALE (L-1) 12b 192 P ya Е E 


To find the marginal density БР 


En Equation (1.74): 


to obtain: 


(1274) 


(X) I integrate over у 


(975) 





L T 
| i 1 
ШИ ж руз |] CL 
LE 4 P (D) 2705 ay L-5 


eo 











2 
207V 


NUES 


where TlT(L) is the Gamma function and K (2) is the modified 


Bessel function of order v. If we make use of the 


identity [20] 


L 
IT - -2+k)! 
K (2) = Eu e - ) eee EN = 
L-> 22 kzl (k-1)! (L-k)! (22)* : 


it is easilv proven that Equation (1.76) is identical to 
Equation (1.67). When x is integrated out in 
Meuation (1.74), I obtain E which confirms the 
result first obtained by Dyer [4]. 

The next step is to solve for the joint density of 
p and p. This is easily accomplished by the following two 


dimensional transformations of Equation (1.74): 


2 0 200 
E. (р,р) = I | | PO. ОАЕ Т) 





go 


which yields: 


er o? be 
P : (р,0) ee > LA ВЕНА 5-7 ІШЕ 78) 
ЦРТЕ ma = E ут 221 2017 


This result is rather remarkable. One can see immediately 
that p and 0 аге independent as they are for the single 
source ES and furthermore that 6 is independent of L, 
and AL the density for L equal receptions is in fact 
identical to FATE the density for the single source. 


Integrating over 6 and p respectively in Equation (1.78), 


P obtain: 





621-1 о? 
P (p £9 СО ЦТ (1379) 
PLE em: е 202 


The mean and variance are respectively 


угп а 


u = 
PLE 2b 


1 еу! 


(L-1)! 
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Mere (21-1) 1! = 1:3.5•,• (21-1) 
2 2 2 
С SL - ми 
PLE 1 PLE 
апа 
Е 1 pS 
| sa (р) = mo eD > (12:50) 
LE 270202 2025 
T 1 -. 
u = 0 and o? = о?у? 
PLE PLE E 


TO complete the statistics for this case, I make the 
final transformation A = e ln x and Á = ex/x in 


Ecuation (1.74): 


e 2 ы 
- (A,A) = sii | | pos ANN (1.81) 
Ap | JAJA LE 


Бо obtaln: 





= О] = 


EA 
| p expl (L*5) 7] "e 12 
ПИНК. |) = EXp/| -e prc een 
Er url ER (n-i): Елі азы 
: (15582) 


Integrating over A in Bauacien (162) торап 


РА (Л) = отт — ехр EL = == exp (2) ШЕ oo 
ҺЕ 2 ау (1-1)! 1 
= 2 = 
DE clin 2g} - y * $102] 


and 


where 





SONS 


til 
О 


L-1 1 
S, (L) = d cee? 5 (1) 
w-l 


| Адаіп this result was previously obtained by Dyer [4]. 


Myeecrating over A in Equation (1.82), I obtain 


. er 1/2v 
P: (A) = == —%; (1.84) 
ALE A a EE 


The moments of A for the equal reception case are 


interesting: 
0 Еп оаа 
: men Ur, 21) 1I (2L-n-2) 1! 
BEAP] = пу 2 "дарын L-n-2)1! | В L (1.85) 
n 
oo ; pou 


Note that the single source case yields an infinite second 
moment as reported in Section l.l, however, for L > 2 


the variance always exists. 
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12.2 Multiple Components with Different 
2115 and Different v?'s 


The most general case of this problem is Case (c) of 
Section 1.2.1, i.e., there are many groups N, and within 
each group there are L, receptions with identical сї апа 
v^. Because there are now two parameters which can vary 
for each source, this implies that in general for 
multiple sources the number of groups will be larger for 
the rate variables than the amplitude variables. I 
define 01. 


receiver pair. Considering the most general case and 


and vi to be oF and v* for the jth source 


Being Equation (1.12) I find the characteristic function 





N МІ: 
О EI 1 
М ШОО = 1. (== (---) ————————— (1.86) 
yX’ X Je eer 71j (w-ilger-+ 202,v202]) 3 
15 ША 
Taking the inverse transform 
| N NL, us ішу ісу 
Bios uud || nn, en 
N е а сы е 1) 


Unfortunately, I have been unable to evaluate this integral; 


thus a retreat for the moment from the most general case 
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is in order. Much more progress can be made in solving 
for the density of X because here we are dealing with a 
sum of independent random variables, and it is not 
required that we know the joint density. This analysis 
will ultimately lead to a very simple and useful 
approximation of the density for x under the most 

general case. Upon completion of this analysis I will 
return to the problem of the joint densities and the pdf's 


of 6 and A. 


1.2.2.2.1 Solutions for Р. (5) 


First I will consider the exact solution when the 


2 
1. 


Base (b).. I apply Equation (1.63) for L= 2: 


product cív is different for each source/receiver pair, 


DEL OS 527-517 
Р. (x) = = = 1 * + e 2 (1.88) 
Хэр 
where b, = 1/20] ,v;v and the * denotes convolution. 
Performing the convolution 
D. bD | EOM 
Р. (у) = Б m poc l e 5 (1.89) 








Tocs 


Rewriting Equation (1.89) in a slightly different way, 


4205 ы bixi by JP -Біхі 
ES (X) - NT 2 телі e = on = e (1.90) 
2р 52-51 b5- bi 
| 22. Әз ТРЗ\Х| | 
Now Р. (у) = Р. (у) % me . Examination of 
Азр Хор 
Equations (1.88 - 1.90) reveals that successive 


convolutions can be done by inspection: 


à M PET 
2 55 1 515, 5, |Х | ba |x 
р. (0) = —— |3 ——— [bye - bye 
X3D OS с 
bi |і 525; -b3 |x] -b3 |x 
ПИ | — > be - b,e (EE OE) 
p2- p? P2 
РЕЛ a 


Rewriting this result so that the method of inspection can 


be applied again, 


x PEE 2.2. ze 


Баба E -b, |x| 
D | 
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bib b, -bjlxl 
DUO. 7 
Bab b, -bjlxl 
Be |_-- 3 (1.92) 


99-7 1-2 2) 4-2 
(Ы; b1) (b$ - 55) 


End thus, 


25354 b; -bilki 
ex x = 2 2 2 2 2 zu 
4D (55 - Бї) (55 - by) (Ба - Бу) 
22272 E : 
E ые: Ж bzix 
2 
D ro ENERO ЕЕ 
(b5 - ру) (65 - 22) (b, - b2) 
223,2 _ 2 
A E oe b3lxl 
2 2 B u 2 
(b$ - bi) (b$ - b2) (bj - b3) 
Bbb: ES 
p 2 2 : | | 2 2 == ) ee 
ЛЕО. . с. 
Thus it can be seen that P. (x) is a weighted sum of 


LD 
the individual single source pdf's of X. It is also 





са 8 


И dent from the L= 2-4 cases that in general 








І 2 
TAID: : 
L ._ j -b. Ў | 
PO (60 =5 ) a u. t (1.94) 
XID | ісі 71 2 2 
mb- = Б) 
m ү 
m-1l 
mzi 
or more simply in terms of 01. and v; 
i о У : 
E. (X) = 2 ) a ехр |- к (05795) 
> ar 1 (or Um go 102) en | 
ТЕТ Та: Im m 
mzi 
Solving for the moments, 
0 | k odd 
- oO L 
L . 2 
k! D dl Iul 7 ЕЕЕ k even 
б i bi L b k+l 
по (Ба - 51) "i 
т= 1 
mzi 


Making use of the properties of sums of independent 
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random variables, 
L 
02: =8 } gv: (1.97) 


Bauation (1.96) with k=2 confirms this result. 
An exact solution for E (X) for the most general 
N 
case considered in the beginning of this section would be 


an N-fold convolution of Equation (1.67) or (1.76) with 


"itself. An exact solution for N=2 will now be presented. 


| Бі (L, +k-2) ! L,-kel 
Р) = Е >, 
2Х p lo kel (k-1)!(L, - k)!2 
L 
2 (L.4m-2) ! I 
1 2 2 
e T ce 95 
(L, - 1)!2 2 m=l Е 2 
“ =b, |y] -B,|x-y! L,-k L,-m 
| о =: ур [ур ау (1.98) 


= CO 


I will now evaluate the integral. I consider ү > 0 


(we know the result for x< 0 must be identical), y» 0, and 








-99. 
then y < 0: 
E» 0 апа у> о, 


Pet y = xt 


L,+L -k=m+l j|" -b,&t -Ь„ў|1-&] L,-k ЈЕ 
| = NS | ЕС en |1-Е | 2 Qt 


О О 


299) 


TO remove the absolute values this intecral must be further 


шта “сеа for t« 1, and t» 1l: 


ECL eLO-k-mel | -bo$ 7! -(b,-5,) Xt L,-x L.-m 
| Am S5 - 2 | EL E лава 
О О 
Mersey -(Ь,+Ь„)ўЕ L,-k L.-m 
Не 2 | oe (1.100) 
1 


The first integral in Equation (1.100) is the confluent 


BUSergeometric function encountered in Section 1.2.1.1. 
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The second integral is a degenerate hypergeometric 
function known as Whittaker's function [18,24]. Performing 
a similar analysis for y« 0 also yields Whittaker's 
function.  Simplifying to as much an extent as possible, 


T obtain for the final result: 


pee d 51105-1 
1 1 72  (L,+k-2)!(L,+n-2) !b,b, |x| 


— rur. A 
2X Di) 212 k=1 m=1 (k-1) ! (L} -k) ! (m-1) ! (5-2) 127777 


Г(о)Г(6) 


‚| ек, FEES 


М[о,0+6, (b,- b.) [Xi] 


од 


expl- 2 (b, -b) |xllte, **pIxl] 2 г) 


+ 


Mas 106 10, +0, 1X1) + (-1)*75 ехр[- 5 (b, - b [Xd d 
2 2 





(048 ) 
e [ке 2 aA x 
E EG. y Ba (dec pO L) 


o — 


à 2 











== 


where, 
а = LQ-ktl 
ô = L,-m+1 


ШІ 0,2) 15 the confluent hypergeometric function, 


and D 12) ssh ttakxer's function. 
, 


Bucrder to check Equation (1.101), if I let 


L4 = L, = Кел I should obtain Equation (1.89). For 
L4 = L, = 1 
„Baba kl 
РО exp(-b,|x|)M[1,2, (b, - b) [x11 
2Х 4 
1 : «14,71 . 
E S b b) ХН Бе БА Wo 2S EO, + ba 1x1] 
1 ‘ ает ; 
* expl- 5 (b5- b) [Xl1 E54 * 55) IXI1 ^ Wo E04 * 55) ITI 


OL 013) 
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лом make use of the following identities [18,24]: 
e? 
Pal, 2,22) = =z sinh Z ! 


Mr UA 


- /z/" = 


where K (z) is the modified Bessel function of order v. 


= obtain 
Wy (2) = М (Z) = expl- 3) 
IA т 002 


and applying these results to Equation (1.10la): 


• 1 е е 
pt = | —— texpt-b, |x]) - exp (-b, 1X1] 
25 em D. i 2 


a [exp(-b, |x|) + exp(-b,|x|)] 





== 


It is easily shown that this result is identical to 
Equation weasel However, though Equation (1.101) may be 
exact, it is of significant complexity and thus of 
limited engineering value. For arbitrary N it is 
therefore clear that Р 0) is analytically out of hand 


at this time. The complexity of the result, however, 


once again motivates an approximate solution. An 


| Edgeworth's approximation turns out to be quite simple 


"— 


and very useful for even the most general case. Recall 


the Edgeworth's series is given by Equation (1.42). In 


the present case, however, because X is symmetric about 
the origin, the coefficient of skew, Үс! is identically 0. 
Thus it remains only for me to solve for ШЕ. the coeffi- 


Event of excess. First I solve for Кі; Ene fourch order 


cumulant or semi-invariant [4 of x: 


2 
ED SEL) (1.102) 


Performing the required calculations on PO given by 


Méquation (1.4) I obtain 


X, - 1920°v* Ci. 103) 


Becalling Equation (1.49) and considering N groups of L; 





ud 


identical sources each, the most general case, I have: 


N 


NEEDS» S Lo... (1.104) 
nó 121 1 ti 


Brom Equation (1.70) I have for N groups of Li Scurces, 


N 
EE N 3L o? v? (1.105) 


 ЕІуіпа Equations (1.44), (1.104), and (1.105), I obtain 


the result 


S S ү,» 

ib E i 
2а 3 I 7 2.1.00) 

( О 0С) 

scu d Bea 

and the Edgeworth's series is then: 
PG) 2 (zn + дут? 2 (1.107) 
X p ге 
М NÁ 


Pr 
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Dee (1.108) 


where п = x/o > and 
NX 


E - exp (- 5 m) 


d. 
von 


It is clear from Equation (1.106) that as N becomes large 
or the L, become large, о approaches a Gaussian. 

The maximum value of Be is 3 for a single source/receiver 
pair and for L identical ones Ya goes to zero as IA NU 
Figure 9 I have plotted the exact density for N=1, [=] 
given by Equation (1.4) with its Edgeworth's approximation 
given by Equation (1.108). This is clearly the worst 
case. I have let gv = l for convenience. In order to 
gain some understanding of the behavior of Vs fOr ENIS 
case, in Table II I have listed rs and Р. (0) from the 
exact density, and its Edgeworth's approximation for 
various values of L identical sources. The percent error 
is also tabulated. From Figure 9 it is evident that the 


error shown in Table II will be the maximum error of the 
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Fig. 9 Comparison of Edgeworths approximation for X when 
N =L = | given by equ. (1.108) to the exact density for 
X given by equ. (1.4). This is the worst case. 
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TABLE II 


The maximum point by point error of the Edge- 
worth's approximation to P. (x) for various 
values of L. ÄLE 

Also listed is the value of Yy., the coefficient 
of excess for each L. S 


e 
Edgeworth Exact ESE 





„2500 22 
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„0938 
s073 L 
·- 0684 
ба 
‚0464 


N л > (09 N HH 
ә UW Б rm DD Ul 
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approximation in the main lobe. From the table 
convergence to the Edgeworth's is very rapid and it 
appears quite reasonable to assume that the Edgeworth's 
approximation will be a very good one when "E $ 1l. It is 
true that the Gaussian assumption will also be reasonable 
For Е < 1, however, because the corrections required by 
Equation (1.108) are trivial, accuracy need not be 


sacrificed for expediency. 


ШИ 2 2 Solutions for the Joint pdf's, 
РО апа Р, (4) 


In this section I will first solve the problem of 
two different source/receiver pairs, i.e., the oF and v° 
are not the same. I will then generale the analysis 
and solve for N=2 of Case (c). Finally, I will indicate 
the analysis required for arbitrary N. 

For two independent pairs, the joint density of 
XırX1rX9 and Xə is given by taking the product of 


Eauation (1.9) with itself with different 01 and v? 


1 
OQ X X2X 7 








Box ox ED RA 
En 2 2 214.19 s 2T4X,0 T 
X 5 X 
es И Ne 1 2 = (1.109) 





LI 


Now, the pdf we are after is the. joint density of x and 

|, X where x = X1 + X2 and y = X1 + Х2- леве та ишла 
 масісп (1.109) X1 and X2 are jointly Gaussian. Because 
| the density of their sum is given by convolution, the 


following can be done by inspection: 


| 1 
р EA 
EE Sx 172 е 
2 402 a Ver (ax, 0777] 1 + 4Х2212У 2 
X X ,2 
К al- E E a | (1.110) 
2077 20% 2(4Х1011У) t 4Х201; v5) 


ИСС final result is obtained using the convolution again: 


X 
u cc A 
Ут [ах 0. > m + &(у-у. DI al 


x 
2D 2005 о 011912 





22 
2011 291 81151191 % Qc Х1)912%5] 


I now let Ху = ХЕ: 





E 


X ?exp(- x/ 207 2) 1 
(0 — —— m 


EB. 
Х›Хәр 807 0? ,v2T , 


2 
2 


LO - c2.) 
= = аа CL ele) 


A 


2995 
2011012 8x0, „V 2 + 8xt (a? N SAU 


I have not been able to evaluate the remaining integral in 
Equation (1.112), however, as will be demonstrated below, 
numerical integration is very simple. Making the 


transformation to p space I obtain: 


Bene, ос 1 
Ер Ò Ep = 2 Jon vo? V5 De ED 2M 
2D 20191972 0 vi + с(011У 1 01295) 
2-52 2 
pug. --0:.) „2 
. 059 «EL, - 2 р Es e TT) 
ES 90 201,9 % 25(611У1-51,У)! 


Integrating Equation (1.113) over p I obtain the marginal 


Mensity of p for two different sources: 





zi 


2 2 | 
ЖИЕ 11212 | 

Р. (0) y 2 2 2 Гао сан та а 52у 

і ? ~ = c 

р E A 02 0 1202 Е о 2 61252) 

с 
E No | Parce 
20127 + 25(511У1 - 01272) 


The integral in Equation (1.114) can be evaluated by 
expanding the exponential in its power series and 


integrating term by term. After significant labor 


E obtain: 


lo о мо) v2 - 02 92) e 


: 2 
P (0) үт (с? P E ы улы l Ср (+) 
2D гуп (94191 - 01,7) n=0 
22 22222 п | 
oe ne? MESES MENS 
ҮНІН г а NEXT д 
- а 7 
ІШІ 25127 ar 2:127 
ESAS») 

where T(a,z) is the incomplete gamma function [23,24]. The 
appearance of || does not mean that P. (0) = 0 because 


2D 
the limit as f>0 of the incomplete gamma function is 


infinite. If the incomplete gamma function is expanded in 





UMS 


its power series, I get: 


(бА - 91 202) со 

(00 89912011 2 
Py (6) = 1 P Cub en) 
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En obtaining numerical results from Equation (1.116) I 
encountered overflow problems in computing successive 
terms of the series before an accurate result could be 
плеса. In order to avoid this problem, Equation (1.116) 


should be rewritten as follows: 
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Values for Br were computed by numerical integration 
пса јоп (1.114) and from Equation (1.117). Depending 
upon the value of ó desired, the numerical integration 

of Equation (1.114) to three place accuracy was about 

40 times faster than use of Equation (1.117). The value 
for pb = 0 is obtained from Equation (1.114) analytically 


and exactly to aid in these comparisons. Performing the 


integration and simplifying as much as possible: 
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The mean of fp is zero and the variance of 6 is obtained 
by interchanging tne order of integration in 


E Eion (1.114), integrating over D first, then t: 
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Despite the complexity of Equations (1.114) and (1.117) 


it appears as though a Gaussian approximation is an 


2 2 2 
ой 
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v5 taken from data, which will be discussed in more 


exceptionally good one. For values of © and 


detail in Chapter III, the values of i from 
Equations (1.114) and (1.117) deviate from the Gaussian, 
using the value of ao Pauar one ел и а Бег 
the fourth significant figure. A more formal 
calculation with equivalent results is obtained if the 
coefficient of excess, м is solved for. I apply the 
Befinition of ME from Equation (1.44 ), perform the 


two integrals required, and obtain e for 9 tor two 


Bifferent sources. 








= з= 


JE bent 6 4 2 
ee, [30_.V, - 30,,v, + 6(0 
è 01, 91 1272 LSI 1 


2 DNA 292 4 
ОЕ и) о с = 365 
|: 1277 Pon Pon 


ale 20 


Applying the values of the parameters used above and 

МЕША Сјоп (1.119) in Equation (1.120) I obtain for this 
example that nr 2.4 x D which certainly warrants the 
Gaussian assumption! Coupled with the result of 

won 1.2.2.1 in which it is proven that Mein remains 
identically Gaussian and independent of the number of 
sources, it appears quite reasonable to state that even 
for the most general cases, Р, (0) will for all practical 
purposes be Gaussian. 


Making the transformation Л = elnx in Equation (1.112), 


hm obtain 
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Integrating over A I obtain 


1 2272 200% 272 а 
не 
р 12[01,v5 et t(o11Vv1 91552] dt 
BÀ E" 42 g? 1 il 1 Su ue qi 5292 ws 22 
| 0 - er tt Bn qM 12"2 4t( E 1272)! = 


es) 


Making an enlightened change of variables allows the final 


integration to be performed, and 
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Though arithmatically messy obtaining the exact solution 
from Equation (1.123) is approximately 50 times faster 
than numerical integration of Equation (1.122) which 
results in three place accuracy. The mean of A is of 
EOurse zero and integrating Equation (1.121) first over 


A and then t, I obtain the variance, 





au 


02 1 6 022 D MENS 

ENERO 17207222112 01201172 5112127 2 312 
Л С? 2 2 4 2 2 4 c 
2D 011912 01» 911 91 2011912 + 211 11 


e күт?) 


Unlike the single source case, the variance is finite 
though the next nonzero moment E[AS.] а5 Infinite. 
The rule derived for equal source/receiver pairs appears 
to apply to unequal ones as well. Namely; ЕЛ5] exists 
only for n even and n « 2L, where L is the number of 


Sources. An example of РА (A) is plotted in Figure 10 
2D 
11! Vie 015: v taken from data. The 


Gaussian using the variance given by Equation (1.124) 


using values for с 


is also plotted. The Gaussian assumption is clearly not 
warranted in this case, and in fact because densities 
for À for all cases are related to the Longuet-Higgins 
type density, the Gaussian assumption will be a poor one 
in general even when the densities are appropriately 
normalized (see Figure 1, Reference [9]), in which case 
the variance of the Gaussian will not be relevant to any 
of the source parameters. 

I will now generalize the approach used to solve the 
two source/receiver case to solve the M=2 case. That is 


two groups of pairs, with Б, pairs each, with the same 
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Fig. lO The pdf of À for two different sources, equ. (1.123) 


is plotted with the Gaussian using the variance 
given by equ. (1.124). 
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01, апа у; (1=1,2). I start by considering the joint 


density of Xy, E, в Xr, E and XL, E given by the 
Beamer of Equation (1.74) with itself with different L., 


01, and vi Eu reni ting it slightly to make the 


Eirst operation obvious: 
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As before, I seek the joint density of y, and Y where 


Х = у + y , and Y = y + ў ~ A double convolution 
ESE LE LE LjE 
Will yield the desired density and because yx and x 
LE LjE 


are jointly Gaussian as before the first convolution 


can be done by inspection: 
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Setting up the final convolution after making a change of 


variables leaves the following result: 
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This is a Key result and is essentially the solution 
of the entire N=2 problem because this joint pdf 
contains all the information required for the solution 
E the pdf's of x, p, A, X, Ò, A. Asa quick check, 
MeL integrate Equation (1.126) over x (which can be 
done by inspection), I obtain Equation (1.20) and hence 
all the amplitude densities for N=2 follow directly. 


Entegrating over x in Equation (1.126), 
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where \ = L,+L_-= 


The remaining integral is clearly not a straight- 








forward one, however, we have already solved for 


P "229 sche previous section given by Equation (1.101) 
2 
and its complexity has been noted. Clearly the X 


densities were handled by much more efficient means in 
the previous section. If I transform Equation (1.126) 
Ento p,f aná A,A Space, than I can obtain expressions for 
the densities P 5 (0) and P (A) and their moments which 


2 2^ 
have been unattainable until now. 
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Entegrating over o I obtain: 
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as expected. Numerical integration of Equation (1.129) is 


we difficult. Solving for the variance, I obtain 
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ШЕ орсатп the final rate density for this case, I 


Meranstform to log variables: 
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and integrating once more over A, I obtain the desired 


| result: 
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Equations (1.131) and (1.132) reduce to Equations (1.121) 
Ed (1.122) for L, =1,=*1. The integral for the variance 
Br f is not straightforward and cannot be performed 
unless Li and L, are known. This completes the analysis 
for the rate densities for Case (c) when N-2. 

The solutions for arbitrary N are very complicated. 
Fortunately, an Edgeworth's series approximation has been 
derived for P g (X) (Equation 1.108) and it has been shown 
that P p ES be accurately approximated by a Gaussian. 
OON MU cannot be successfully approximated 
by either an ED cd. S Series or the Gaussian. The 
coefficients of skew and excess cannot be found for À 


Bor arbitrary N, thus ruling out the Edgeworth's series, 


and though one can certainly fit a Gaussian to NES 
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шшосеа earlier, the fit will not be good, and the variance 
_ obtained from the Gaussian will not be relevant to any 
physical parameters of the problem. For realizations 

with N>2 in which P NO is needed, and/or approximations 
For Es Or Py) are not good enough, the following 
procedure can be used for arbitrary N provided computer 
time is available. 


For arbitrary N, following the approach used for 


| Ш- 2, І сап write: 
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As before, the sum of the XQ S involves the convolution of 





ee oe 


Gaussians which can be done by inspection and 
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Unfortunately, this is as far as the analysis can be 
taken without considering a specific N. Thus, the 
procedure would be to numerically compute N convolutions 
with respect to the X,'S to obtain PARA Then one 
more integration over x yields the exact density for 

Xy likewise, integration over p and A after simple 


transformation yields the exact density for бы апа Ax 


respectively. Perhaps an easier route to P И ~ (Хи) 
N / 
can be taken by numerically performing the two 


integrations required by Equation (1.87). 
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In summary, the solutions for the rate densities 
for Case (c) with N>2 are of sufficient analytical' 
complexity that approximate methods or numerical 
techniques must be employed. The densities for even the 
most complicated realizations of Case (c) for 6 and x 
are well in hand with the Gaussian and Edgeworth's 
series approximations respectively. For A, however, 


complicated multiple source configurations must be 


handled numerically. 


e’ Coherent Source Addition 

In this section of the thesis, I investigate the 
statistics of the received signal when two or more sources 
are radiating at the same frequency or so near in frequency 
that the receiver is unable to resolve them. The 
application of this analysis is therefore dependent upon 
parameters of the receiver; i.e., its resolution R, which 


is a function of its averaging time T (R YT 1/T), and 


the source bandwidths, specifically R> Af, and R?Bi,Vi, 


| 





Bere Af is the frequency separation, if any, of the 


Sources and В, is the source bandwidth of the ith source. 
With these conditions, the procedures used in previous 
sections are not applicable, and coherent addition 


of the signals must be considered. This analysis will 
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model random, narrowband, sea noise from biologics, 
weather, distant shipping, etc. that overlap in frequency 
in a small band of interest. 

The multipath signal, p(t), under the assumptions 
of phase random propagation can be written for a single 


source a5 


N 
БИО = = ) Соз (ш + Sad 2135) 
n=1 


where r is the single path amplitude, N is the number 
of paths, and 8, is the single path phase, which is 
distributed uniformly between 0 and 27. For many 


sources at or very near in frequency as indicated 


above, 
È N. 
с i 
КОО = )| г. ) cos(ut * 8 ) (1.136) 
1 n 
1=1 п=] 
where Lo is the number of coherent sources. If this 


analysis is applied to random sea noise, then 


Bouations (1.135) and (1.136) imply that each radiator 








is CW with essentially constant mean output levels. This 


is perhaps not too bad an assumption if integration 


times and record lengths are not very long. It is 

assumed that N: 24, Vi. Such that phase random propagation 
is obtained for each source. Thus, Equation (1.136) 
reveals that the problem is one of solving for the 
Statistics of a random vector which is the sum of Lo 
random vectors, each with Rayleigh distributed magnitudes 
with different means, and uniformly đistributed phases. 
Forming the envelope of Equation (1.135), I obtain the 
quadrature components for one source, and applying 
Equation (1.136), I obtain for the quadrature components 


of the total vector: 
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where p is the rms amplitude of the total vector, ó is the 


total multipath phase, Ds is the rms amplitude of the nth 


Source, and $a the multipath phase of the nth source. 
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The quadrature components for the single source are 
Gaussian [8,9], thus from Eqn (1.137) it is clear that X 


and Y are Gaussian as well: 
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Бе that the short hand notation for Y and Y is meant to 


Show that the form of the densities is identical. The 


quadrature components are in fact independent. Eqns (1.138) 


and (1.139) illustrate that no conditions need be placed on 


the strengths of the sources or relative values of the 


Ја '5. At this point, however, I should caution that if 


One or two sources are propagating energy via 3 or fewer 
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paths, and they are strong relative to the other sources, 
dominant path affects must be considered. Please refer 
to Appendix C, and for dominant paths plus diffuse 
(Rayleigh) noise see Reference [26]. 

Continuing from Equations (1.138) and (1.139), the 
analysis is straightforward. The resulting densities are 
meet identical in form to the densities for the single 
narrowband source as given in Appendix A. The functional 
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| difference is that wherever c2 and v? appear in Equations 








Т 
(А1-А8) they must be replaced bv oF and 011/95 respectively, 
B.q., 
5-9 Y 
| Р (О) = О expl- 0” > 0 
| біз О Ел o p 256 , О 
с 
қы v2T 2 i 2 ТІ 
E cm с. : 7 = 01(2- 5) CESSIT) 
L L 
с с 
ПИ : ___-- 1 у. 
Ar de [1 + —=— ( —L.) 121 3/2 { 
С ác? Orr 
Ш- = 0 Р of = 2 8 141) 
| Ar Az, 


| С С 
| 
I 


| The most important difference, however, is that now the 


densities for Л апа # depend upon the o? 


in 


| Single source they do not. It should be noted, however, 


's while for a 


Menat if v^ is the same for all the source/receiver pairs, 
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then once again A and p are independent of the amplitude 
parameters, 011158; and are in fact equivalent to the 
densities of A and Ф for a single source. Even though all 
the sources are at the same or near in frequency the 
likelihood that v“ is the same or almost the same may not 
be a good one given v*'s range dependence [11] (when the 
internal wave model applies, see Section 3.2). 

I now consider a narrowband experiment in which random, 
background sea noise is present. From the analysis above, 


I can write: 


and 142) 
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р 1 YT 
п=1 
where Lo is now the number of noise sources in the analysis 
band and 01 ӘП о are the signal parameters, and the din 


and у are the noise related parameters. I now make the 


assumption (perhaps bad) that v? is a constant for the 


Signal and the noise. Thus, E = о; and I define 
Lc 
e S nr which is one half the sta mean square noise, 


O 
and therefore: 
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where SNR = Ong time average mean -guare signal | 


long time average mean square noise 
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Applying Equations (1.144) to the densities and statistics, 

e.g., Eqns. (1.140) and (1.141), the effect of background noise 
ene accounted for. It is clear that SNR will be a function of 
the analysis bandwidth, and that noise levels should be 

measured in the absence of the signal. By narrowband I mean 

the signal bandwidth B should be <<2v which in turn should be 

<< 1/T where T is the "E time and the analysis band- 
width ВА should not be so large that the frequency dependence 


BF v 


Would cause v^ to vary significantly across the band. 
If for other reasons, i.e., range dependence, the 
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Finally, if one is analyzing over a wide band and the 
components are disjoint but each is a sum of one or more 
sources at that frequency, then all the results of the 


frequency disjoint analysis can be applied with 91 replaced 


2 
I 


defined as given in Equations (1.138) and (1.139). 


By c2 and v?^ replaced by roe where these parameters are 
For noise which is continuous over the analysis band, 
the statistics will be a function of the total bandwidth 
and the receiver resolution. All the analysis of coherent 
and incoherent (or more correctly, frequency disjoint and 
independent) sources apply except now the number of groups 


N is given by B,/R where B, is the analysis bandwidth and 


A 


R is the receiver resolution. Because the sources within 


| and en must be found for each 


ons using Equations (1.138) and (1.139). The groups can 


each group are coherent c 


now be considered for purposes of applying the results of 
Sections 1.2.1 and 1.2.2 as individual source/receiver pairs, 


Eth the o2's given by the o2's and v*'s given by the 


1 I 
E/or's. In conclusion, therefore, the statistics of the 
received signal are a function of receiver parameters as 
well as source parameters. For large observation or analysis 


| Біте Т, and consequently higher receiver resolution, the 


EncGlysis of Sections 1.2.1 and 1.2.2 will most probably apply. 
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As T becomes smaller, however, coherent effects must be 
considered as developed in this section. Careful regard for 
the inequalities stated in the beginning of this section 
must be applied, case by case, to the data set of interest 
and the receiver characteristics, in determining which 


analysis applies. 


m2. 4 Crossing Rate Statistics 


In previous sections I have solved for the joint 
densities of amplitude, and amplitude rate for many multiple 
source cases. This 
allows me, therefore, to solve for the theoretical 
amplitude crosSing rates for these cases as well. 

Following Rice [10] and Dyer and Shepard [15], the 
mean crossing rate for the sta rms pressure о сап be 


defined as 


со 
С(р.) = | |бір, (0,0249 (1.146) 
= со 
where Po is the axis crossing level. For L equal 
source/receiver pairs, Equation (1.78) is used in 


Ecuation (1.146) and I obtain: 
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Bauation (9) of Reference [15] is obtained for L=l. 

(Note that р. = оѓ, р. = 6292.) For two different 
О ДЕ 2 1 

sources/receptions Equation (1.146) is applied to 


EStion (1.113): 
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If I make the change of variables 2 = a+öt the integral in 


Equation (1.148) can be evaluated: 
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where y(v,u) is the incomplete gamma function [23,24]. 
Finally, for Case (c) N22 I apply Equation (1.146) to 


Equation (1.128) and 
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where a, 6, and c are the same as above. 


Integrating I obtain [24]: 
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where T(x) is the gamma function [24], and $1 (0,8,Y:x,y) is 

a degenerate hypergeometric series in two variables [24]. 
For coherent sources, phase crossing rates as well 

as amplitude crossing rates can be found. As with 

amplitude, for the multipath phase [10,15] 


со 


с(%2) = | |ФІР, 29,94 (1.152) 


= со 


Following Dyer and Shepard [15] but allowing for the 


effect of many sources at the same frequency, it is easily 








shown that 
ОТТ 06 
сы (о) = y2/T ce expí- p m 153) 
С T I 
and 
с 
ee re: 
бү Te PEE ОТ (6219520) 
e 
Thus, 
G(p_) „ О ја 
ec 4 "Oo O 
CT.) = (87) о exp ( 205 6,155) 


where 01 апа m are given below Equations (1.138) and 


(1.139) respectively. Allowing for the normalization 





EI 


2 2 
I and OTI! 


Bue form of Equations (1.153 - 1.155) are identical to 


employed in Reference [15], and the parameters c 


Equations (9) and (10) in Reference [15]. It is also 


noteworthy that Equation (1.153) remains 


21 


unchanged whether the v“'s for all the coherent sources 


are the same or not. In terms of the SNR then 





mo (am)? (SNR_,[0 e % ENS (1 156) 
Ge) пена со BLZ 75. NT SNR 


I now consider the ratio of Equation (1.155) with 
Equation (10) of Reference [15] (the limit of Equation 1.156 


as SNR > о), 


2 
R= 5 solz T (1.157) 
Thus the effect of noise on the ratio of amplitude 
crossing rates to phase crossing rates is a multiplicative 
factor given by Equation (1.157) which depends on both the 
SNR and the axis crossing level selected. 
Finally, if one desires the crossing rates for the 


mean square pressure y or the level in decibels A, the 


О n“ 


simple substitutions C and ade ехр(А /2є), respectively 


is all that is required. 





-142- 


1.2.5 Summary of Results 

Because of the many cases, the notational difficulties, 
the complexity of the equations, and the plethora of 
random variables involved, I will attempt a summary in 
the form of a table of the multiple source section to aid 
the reader in gaining a little perspective on what exactly 
has been accomplished. In Table III, I have compiled 
the overall results of the analysis for the various pdf's 
of interest in the multicomponent and single source cases. 
On the left is the breakdown into cases based on source 
structure and signal analysis. Across the top are the 
random variables whose pdf's we seek. When a number 
appears alone in a box, it indicates the equation number 
Of the pdf and also that the result is unique to this 
thesis.  Superscripts (circled numbers) appear in boxes 
for comments below, and numbers in brackets refer the 
reader to those references in which the equation appeared 


previous to this work. 


1.3 Finite Bandwidth and/or Modulated Source 

In many actual oceanic situations the source does not 
exhibit stability in frequency but in fact oscillates or 
wobbles about a center frequency which can be characterized 


by either frequency or phase modulation. The effect of 
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this modulation on the family of single source densities 
is analyzed. Even "narrowband" signals have finite though 
small bandwidths. I have determined a criterion for 
smallness, as well as the bandwidth effect when this 
criterion is not met. Furthermore, the analysis presented 
in this section reveals a method by which finite 
bandwidths and/or source induced modulation can be 


determined from the received signal. 


1.3.1 Amplitude and Amplitude Rate Densities 
In the absence of modulation we can write the multi- 
path signal, p(t), under the assumptions of the phase 
random model as given by Equation (1.135). When the 
source is frequency or phase modulated, we can write 
p(t) as 
N 
EU  cos[ot - M(t) - 9.] mus 
п=] 
where M(t) is a function of time which may be random that 
represents the modulation. As indicated in Equation 
(1.158), it is assumed that source induced modulation 
will be path independent which implies that any change in 
signal propagation characteristics (i.e., path structure, 


or volumetric absorbtion) will be independent of the 
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"instantaneous" carrier. This is apt to be so unless the 
modulation is extreme. This assumption also applies to 
the finite bandwidth effect. The bandwidth cannot be so 
large that different propagation characteristics obtain 
for the extremities of the signal. Further, it is 
assumed that none of the energy in the signal is rejected 
because the "instantaneous" carrier is outside the 
bandwidth of the receiver. Likewise, for a finite 
bandwidth source, it is assumed that the entire signal 
bandwidth is within the bandwidth of the receiver. Note 
that depending upon the specific temporal CSS E 
M(t) the modulation would be classified as either frequency 
modulation or phase modulation. This distinction, 
however, does not alter the analysis to follow. 

I perform quadrature demodulation on p(t) and obtain 


the quadrature components: 


N 
X = pcos[M(t) + 6] = xr ) cos[M(t) + 6] 
n=l 
MR ESO) 
N 
Y - psin[M(t) + ф] = r ) sin[M(t) + IR 
n=l 


where p and ¢ are the amplitude and phase respectively of 
the complex envelope, or alternatively in the terminology 


of phase random acoustic propagation the sta rms pressure 





ao 


and multipath phase. Taking the derivatives with respect 


to time, 
N • Б 
X - -r ) [M(t) + 9 ]sin[M(t) + 9 ] 
п=] ý b 
(10160) 
| М. | 
Y=r) [M(t) + 8 _ ]Jcos{M(t) + 8.] 
n-l à d 


In terms of the quadrature components, Equations (1.159) 


and their derivatives, Equations (1.160), I can write: 


(x2 + y2)% 


© 
|| 


Cr on 


EE 
2 


De 
| 


(XX + YY) (Xx? + y?) 


ШЕ civen in Section 1.1, x, x, A, and A can be expressed 
Ben terms of р апа р. By applying Equations (1.159) and 
Ш 2169) to Equations (1.161) and making use of 
trigonometric identities, it is easily proven that 
Equations (1.161) are independent of M(t) and, in fact, 
are equal to the result obtained when M(t) =M(t) =0. 
Thus, the amplitude, and amplitude rate variables are in 
fact independent of the modulation. It also follows, 
therefore, that the joint densities of amplitude and 


amplitude rate are independent of the modulation. 
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1.3.2 Multipath Phase Rate Densities 
For the multipath phase rate, the modulation plays a 
Bratical role. From Equations (1.159) it is clear that 


the multipath phase with modulation, Фм’ is given by: 


E x 
On DE x^ Meyer % (le 162) 
and, therefore: 
à, = M(t) + 4 (1.163) 


This result can also be obtained from the single path 
variables alone [the extreme right-hand Side of 


Bouations (1.159)] from which, 


N N 
> E r? E . 
9м = A pry |. соз [M(t) + ИЯ М (+) ng 9 12059, 
N М. | 
+ } sin[M(t) +0 ] ) [ice + 5, 1sine, | (1.164) 
п=1 n=1 


Again, using trigonometric identities, I find that 
Equation (1.164) is equivalent to Equation (1.163), in 
which ó is given by Equation (1.164) with M(t) =M(t) =0. 


mote that from Equations (1.159), Equation (1.164) is 








zn 


independent of r. 
M(t) can be either a known deterministic function of 
time or a random process governed by a pdf EE a 
the latter case, because M(t) is independent of ocean 
Mearameters, Equation (1.163) reveals that Фм is the sum 


of two independent random variables. Therefore [14], 


mda) = p. (MY&D- (à 
Pj (0) Pg (M) *P; (0) (1.165) 


where * denotes convolution. The mean u = ру, and as 


ф 
М 
with an unmodulated source, the variance 05 2 
M 
I£ M(t) is a known deterministic function of time, 
fr define: 
Ts 
Я | 1423) 
(0) - =l >> ДЕ [NE 
Фм ДЕ гт + [6 = M(e)]?,:72 
0 y 


where ТЕ is the length of the time series (not to be 
confused with T, the averaging time of the receiver). 


| 


CMM is a continuous histogram and has all the properties 
Of a pdf, i.e., it is always positive and integrates to 
| one. This function or pseudo pdf can be employed when 
| M(t) is deterministic but not periodic for a given ensemble 


of time series. Equation (1.166) also applies for 


periodic deterministic M(t). However, for periodic 





E 


deterministic M(t) puwhrich there is exactly опе ог п 
integral number of periods M(t) can be randomized, treated 
as if it were a random variable, and its pdf found 
enabling use of Equation (1.165). If there are many 
periods in the record, then an integral number is not 
required; however, some error will be introduced. As will 
be demonstrated by examples below and in Section 1.3.3 

for periodic modulation functions, Equations (1.165) and 
(1.166) yield identical results. For many interesting 
problems in the ocean, M(t) may be deterministic but 
unknown, the real (nonrandom) parameter estimation 

problem [25]. In these Situations one will obtain 
experimental realizations of н (9 from which one is 

able to learn characteristics of M(t), as will be 
demonstrated in Chapter 3. Also, as will be shown in the 
figures, some deterministic modulation functions will have 


easily recognizable histograms, Н: OU and in fact 
M 


Ф 
(6) is by itself a valuable piece of 


^u 


Enrormation to have. 


knowledge of H 


I:shall now consider three analytical examples 
illustrating the effect first of sinusoidal phase 
modulation, second of uniform frequency modulation, and 
third of Gaussian frequency modulation on the pdf for 5. 


For sinusoidal phase modulation, 








== 


M(t) = Bsin(ot + > 5) 
and thus, o 157) 
M(t) = Bocos(ot + фы) 


I have randomized the phase witn dua ТЕТІК стер авт сеп 
between 0 and 27), which indicates uncertainty in initial 


conditions. I obtain the pdf for M(t) [3], 


1. 


——— ———— , |M| < Bo (1.168) 
c = 42) 1/2 | 


: p i 
PgUD = т 


Ка плод Equations (1.168) and (A8) in Equation (1.165), 


B find 
2 Bg 
Pa (х) = 57 | A O OA (1.169) 
M ае А ает у, 


An analytical expression has not been found for this 
integral; however, numerical integration is straightforward. 
Applying Equation (1.166) with те 2тп/с (п іѕ апу 
integer) and M(t) as civen by Equation (1.167) with %_ = 0 


(don't forget to exploit the symmetry of the cosine) 








| 


| 


Bb 


yields Equation (1.169) as well. When x=0 in 


Equation (1.169), the integral can be evaluated: 


P (0) = a E с ОНИ 2.170) 
M ту: + 8:02: 1/2, vs: +B°o° 


where E[T/2, k] is the complete elliptic integral of the 
second kind. 

It is possible to make some progress in solving for 
5, 5) if I make use of the convolution property of 


Fourier transforms or characteristic functions. For 


independent random variables [14], 
| = u) ELE) 


where М, (ш) is the Fourier transform or characteristic 
function of the pdf of random variable x [14]. Шү) 15 


given in the Appendix under Equation (A8), and for the 


sinusoidal density Equation (1.168) [24], 


Mg (w) = 22 (980) 


where Jg (2) is the Bessel's function of zeroth order. 








EO 


Thus, by exploiting symmetry and taking the inverse 


Eransform of Equation (1.171), 


со 


ES. = ~ | uR, (wv)J_(wBo)cos шх ах (MA) 


0 


where K, (2) is the modified Bessel function of order one. 


Expanding the cosine and integrating term by term [24], 


со 0 2n 2 
P (x) = = ) — (=) (п + 5) (020-1) 11) 
М п=0 ` 


2n * 3 2n * 1 , Par 
=o , De / 1 , кы Ке 3) 
een - |) :! = 1•3•-5 ... (2n- 1); (-1)!! = 1, апа 


F(a,b;c;z) is Gauss's hypergeometric series [24]. 
Unfortunately, Equation (1.173) converges only for 


Bo<v and x<v. When Bo=v, 





= Por 


P(x) = /177 2 Y ED m Ltt 
Ж = ТІ - = = n- 11 
Фу uo 
i 1 
S d——  — — I MH (1.174) 
noct n DM me 
G a GSS SAA йт 


ms With Equation (1.173), Equation (1.174) converges only 
FOr x< vV. Р; ($) is plotted in Figure 11 for the case СЕ 
no modulation, Equation (A8), and for various values of 
$c relative to v using Equation (1.169). Applying 
Eguations (1.170), (1.173), and (1.174), when applicable, 
revealed that the error of the numerical integration of 
Equation (1.169) is approximately 1%. 

For uniform frequency modulation, M(t) = Mt and 


M(t) =M, where M is a uniformly distributed random 


variable: 


NADA |М| <А 
(М) = (Ty o) 


0 otherwise 


and A is the maximum excursion from the carrier in Hz, 
thus, 2A is the bandwidth of the modulation. This 


Bharacterization of M(t) would apply, for example, to a 
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deterministic M(t) such as a saw tooth. As before, I 
ШЕріу Equation (1.165) using Equation (1.175) and (A8), to 


obtain: 


P^ (x) = а eee (1.176) 


ОМ NE RFA WFAN? 


Euuation (1.176) is plotted in Figure 12 with no modulation 
(Equation A8) and for various values of A relative to v. 


For the final example, Gaussian frequency modulation, 


m2 
| -% 
PM) = —=— е 96 (1.177) 
ү ото, 


where e is the variance of the modulation. As before, 


m convolve Equation (A8) this time with Equation (1.177) 


and 


со 3 
2 => == 
| 2 
Pa (xX) = ~~ | (v? ty?) * exp[- == 
M Стос 





(х- у) 2]ау 


GN 


ОЕЕО) 


As with sinusoidal phase modulation, I have been unable to 


solve the convolution integral analytically. However, 
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о улла the characteristic functions, I obtain an 


alternate expression of Equation (1.178): 


> = У p 
B T | we K4 (ov) cos wx do 179) 


Expanding the cosine, I obtain, 




















oo n 2n 
a - 2n + 2n + 
P w = — } 435 (e CERA 
M с-ту2 n=0 ` G 
E E 
* exp е 20) W 1 ı (377) (1.180) 
Eoo. 
Or, alternatively, expanding the exponential in 
Beuation (1.179): 
ИШ ®1 ; 1 (26927 pata pati 
фм El = ту КО v 2 2 
2 
P (2 үші, si > SS) (le 1 $1) 





where 172) is Whittaker's function, and F(a,b;c;z) is 
Gauss's hypergeometric series. 

шта ју, Equation (1.176) will also model the effect 
of a nonmodulated, but finite bandwidth source if the 
conditions stated in the beginning of this section hold. 
It is clear that as long as the energy is uniformly 
distributed on the average between Е A and Е ФА. 
where fo is the carrier frequency, then Equation (1.176) 
applies, and the bandwidth, B, is given by B=2A. A 
glance at Figure 12 reveals that the effects of the 
bandwidth on the pdf for ф сап ђе neglected if B«« 2v. 


The above analysis also reveals that modulation coupled 


with bandwidth effects are additive. Thus, 


E = ААА (OU 


where Po (ф) is the pdf of the bandwidth which we have 
B 


assumed is uniform between +B/2, and Foam A is the pdf 
For 4 when a narrowband signal is modulated and the 
Criterion B<< 2v is not satisfied. 

Taking the three examples of modulation used above, 


and coupling them with the bandwidth effect, I obtain for 


Sinusoidal phase modulation, 








Bo 
er _ 
MB a 8°с° - y? yv? - (x-B/2-vy)? 
B ШЕЕ y dy 


yv? + (x - B/2 = y)? 


or, 





• = 200 = 1 wB 
| > ELE | K, (wv)J, (w80)cos UX SIN 5 do 
О 
For Gaussian frequency modulation, 
1 I 2 
E (x) = — ехр[- —r (x-y)“] 
Фив 2BY 27104 m 206 
. y * B/2 E v = B/2 |, 
NAFTA era ae 
or 
= - 20262 
P' iss 2v | e K. (wv)coswx sin E du 
PMB BT l 2 
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(1.184) 
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End for uniform frequency modulation, 


Р 20 = == [^n + (XtA+B/2)* - YVv* + (x+A-B/2)* 
+ fy? + re 2—04 у^ + GENS 
· (1.187) 


Finally, it should be noted at this point that because 
a finite bandwidth signal is indistinguishable from a 
"narrowband" signal which is experiencing extremely rapid 
uniform frequency modulation as indicated above, the effect 
of the bandwidth, as with modulation, is felt only by Ф, 
the phase rate, and the amplitude variables remain 


unaffected. 


IM Crossing Rate Statistics 


For path independent source induced modulation, I 

|, have shown that the amplitude, amplitude rate, and joint 
densities of amplitude and amplitude rate are unaffected 
| by the modulation. It follows, therefore, that G (59) 
will be independent of the modulation and will be equiva- 


lent to the nonmodulated single source result given by 





Elo 


Beuation (1.18). 


For a modulated or finite bandwidth signal, ene! 


is as expected critically dependent on the parameters of 
the Pedulation and the bandwidth. To find G ($9) ЕЕЕ 
necessary first to solve for Py 510,9) and this 501 пена 
is crucially dependent upon the exact nature of the 
modulation or bandwidth. For deterministic modulation, 


EEconsider Equations (1.162) and (1.163). Тһе joint pdf 


for $ and ф without modulation is given by [8,9]: 


: 0«6«2m (1.188) 


> = 2 


52 272 
DE, 22! 


Making the change of variables given by Equations (1.162) 


ЕСІ 1.163) I obtain 


is 
m ф ono) = rar | — 375 aie (1:139) 
M S ; [1 + 42 — MN) 


For the interval 0 to то: min[M(t)] «40 « 2r * max[M(t)] 
апа |ф| < ©. Note again that the H function has all the 


properties of a pdf. To find the phase crossing rate 
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Gx (Po)? apply Equation (1.152) and 


т 


со 


| | ње ало 


END 5 LI——-—— ; : 
Mo 4туТ Ж 2/42 


Rewriting Equation (1.190), I obtain 


T 
2 S : 
М Ф Ў 
Де) = {an | | Br ы Т оу. т dé 
E ATT. „102 + M(t)? € 28 (6)6 € 82] 7* 
0 
E | A 37 dé dt (1.191) 
[v? +M(t)? - 2M(t)6 + 6?] 


Performing the integrations over à: 





е) = == | /у2 + M(t)? dt (1.192) 
5 


This is as far as one can proceed without the exact form 
ШЕ M(t). (Note that if M(t) > 0 I recover the no 


modulation result.) Equation (1.192) can also be used to 





AA УУ | ы 


EMO 


find the phase crossing rate when nonstationarities make 
у E function of time as well. Obviously, however, v(t) 
must be known. If M(t) is a random process or the 
modulation is periodic and exactly one or more integral 
number of periods are on the record (if there are many 
periods the number need not be integral), then the 
probabilistic approach used in the previous section can 
be applied here as well. From Equations (1.189), (1.165), 


and (1.166), 
. 1 • . 
P, 5| (06,0) = = Pi (6) (1.193) 


For sinusoidal phase modulation, I use Equation (1.169) 


Bi Equations (1.193) and (1.151) to obtain 


с ( m E J : х 
ks! de] ee | nn MY 
А) 
Integrating first over x I obtain 
| Bo а 
Gulo) zl: | AES =e ye ey к КОЕ) 
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The remaining integral is a complete elliptic integral 


of the second kind: 


СФ) = I. RF EiS , r] ПОО; 
where 
o c 5 _ 
GF? 


For either В ог сос ~ 0, См (%6) for no modulation is 
obtained. 

To demonstrate the equivalence of the probabilistic 
approach and the deterministic approach given bv 


Equation (1.192), I consider again sinusoidal phase 


modulation: 

M(t) = ßsin(cot + ф <) 
апа 

M(t) = ßocos(cot + Фф.) 


Assuming I have exactly n cycles, 





| vv? +8?7°cos-ot dt (00197) 





ee 


L observe that the integral over one cycle of a cos*x is 
equal to four times the integral over 0 to 1/2. Applying 


this observation to Equation (1.197), 


1/20 
_ 1 4no 9 9 — 
Gulo) or элт | INS COS Oe die 
O 
and therefore 
1/20 
Gu ($) = = | НИН oe ge ce (1.198) 


I make the change of variables x= с+ and apply the 


trigonometric identity соѕ?х = 1 - sin?x in 


Equation (1.198) and simplifying, 


1/72 


Gy (to) = = IS ag | vI = [8°с6°/(%° +8°с7)]зїп?х dx 
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The integral in Equation (1.199) is the definitive form 








Er the complete elliptic integral of the second kind and 
Equation (1.196) is recovered exactly. 
For GausSian frequency modulation, Equation (1.178) 


is applied and 





у^ | 1 | 1 
ШЕ) - -- — | x expl- +> (x - y) ]&хӣу 
ү о 21 2106 (v2 + y2) 3/2 206 
E О 
У) 
Performing the integration over x yields 
«уль 
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= G M V 
G(¢.) = ——= er Da) 4... (--) (ду 
М ‘о (27) 3/2 (у2-ү2) 3/2 2 m 2 G 
32015420 1) 
where Е istene parabolic cylinder function [24]. I 
now apply the identity [24] 
T Zunge B 
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p lzp 202272 р Zu? 2 
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zo е 


о Equation (1.193) and I have 
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ю ехр(- a 





5950 2 
G0.) = З Е 1,5,2) бу (1.202) 
1/20 (v? +y?) бе 


where M(a,b,z) is the confluent hypergeometric function. 
The final integral in Equation (1.202) must be performed 
numerically. 

For uniform frequency modulation, direct use of 
Equation (1.176) yields a value for сс Ее. This 
result is not a physical characteristic, but rather a 
И пепсе only of the mathematical form of 
mGuation (1.176). The integral leading to Equation (1.176) 
is 
102203) 


ee а 
2 оа 
-A ату) 


ЕУ Баласіоп (1.203) to Equations (1.193) and (1.152) 


yields 





IU OE 


I make the change of variables t = x-y and 


р! tty 
Eu) = aA | | 
d 


SS dt dy 


-A - 


Performing the integration over t yields 


2 а 2 
V ШІ v 
M o 4TA 2——2 — 2 /—2———2— 
la у у ууу“ + у 


and finally I obtain the result 








ns V Jur ee 
G(¢_) 9 —— ln|————————| + ә ио + А (19204) 
М о 4TA E DE = - ar 


In the limit as A + 0 Equation (1.204) converges to the 
no modulation result. Note if I let A=B/2 in 
Equation (1.204) then I have exactly the phase crossing 
rate for a finite bandwidth non-modulated source. 

The crossing rates for modulated signals with band- 


width have not been solved; however, the procedure is quite 





zo 


straightforward, though the integrals may not be. 
I now consider the ratio of amplitude crossing rates 
to phase crossing rates. Dyer and Shepard [15] obtained 


this ratio for the narrowband, no modulation, single 





source: 
G (p) y PL 0- 
= (817)? == exp (- 322) (1.205) 
С (Фф) 1 1 


Equation (1.205) is independent of v and depends only upon 


2 
i 


is a controllable parameter unrelated to oceanic phenomena. 


ci which, being a measure of the energy in the signal, 
However, non-stationary behavior of 91 due to ocean 

i namics will affect the ratio given by Equation (1.205). 
Likewise, if the source is modulated or has a bandwidth 
which is not << 2v then the ratio given by Equation (1.205) 
Will be affected. For these cases, the ratio G (05) /G (95) 
ШОО be a function of c2, v*, and parameters of the 


modulation. Adopting the approach in Reference [15],. 


I obtain 


MP) 
41%) 





= % 
= (8m) Coi P. (05) 12206) 


where, for sinusoidal phase modulation, 
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ЕЕЕ a Bo B (1.207) 


2/8?05?FV?2 Ely,r] ти 





for uniform frequency modulation, 


с = 02200) -—  _ (1.208) 


T X ] + Ау)? + А 





UEFA? - А 


for Gaussian frequency modulation, 


-1 
БЕ 1525004)! (1.209) 


where G (60) is given by Equation (1.202), and £inally 
for a finite bandwidth, non-modulated source C is given 
by Equation (1.208) with A=B/2. Except for the finite 
bandwidth result, great care should be taken in applying 
the formulas in this section to insure that the actual 
modulation fits the kinds of modulation assumed here. 
However, the procedure developed in this section can be 
applied on a case by case basis to any kind of frequency 
or phase modulation to determine the phase crossing rate 


factor. 
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CHAPTER 2 


COMPUTERSSTMULATION 


A computer simulation of phase random multipath 
propagation was originally developed for two reasons. 
First, in the course of analysis of the multiple source 
cases it was felt that a computer simulation would provide 
confirmation of the rather complicated analysis when 
datawere unavailable or difficult to obtain for the case 
in question. Second, when an PI impasse was 
reached the simulation could provide insight into the 
nature of the solution, thus aiding in the analytical 
process. Тһе simulation fulfilled these two objectives 
not only as originally intended in the area of multiple 
sources but in all aspects of the theoretical development 
presented in this thesis, as well as providing confirmation 
Of some of the basic precepts of the phase random model 


of acoustic propagation. 


2.1 Computer Model of Phase Random Multipath Propagation 


A computer simulation which generates random samples 
from a phase-random multipath process has been developed. 
Currently there are two versions. The first, called 
RANDPHASE, simulates up to 50 equal intensity, equal v°* 


components, including sinusoidal phase modulation, uniform 











-172- 


frequency modulation, and bandwidth effects on phase rate 
for a single source. The second, called BURMRAN, is 
capable of handling an arbitrary number of groups of 
Enecual intensity but equal v^, for the amplitude and 
amplitude rate variables. Both programs are written in 
FORTRAN IV and were run on an Interdata Model 80 computer 
with an IMLAC display processor. The FORTRAN listings 
for RANDPHASE and BURMRAN are contained in Reference [27]. 
The following algorithms are applied to generate 


samples for the amplitude and amplitude rate variables: 


2 2 
allt | 

К = к: cosa -- Sing 

E ie 1-41 Pi =1 1 
ГІ 

К = 285 sind 8 cos9 

1 Lu A a a a 

N N 
- ) cose ) 6 sine 
п=] Пі n=1 1 un 
he 17/2 

E- X 

ei ,-1/72, 

E 2x XL 
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> 
| 


L = 10 10410 Xr, 


> е 
| 


а= [10 log, ge] Хг XL (2.1a-£) 


where г, is the single path amplitude for the ith source, 
L is the number of sources, N is the number of paths, 
9, the single path phase, and Т the single path phase 
Ee. For L=1, the single 2. samples of the 
multipath phase ọ and the multipath phase rate dé are 
generated. Equation (1.164) is used to generate the 
samples of $. The input parameters to the program are 
eeu, N, A (or equivalently, B/2), 8, c, and I, the 
number of samples desired. For RANDPHASE the r, are 
equal and are set to one for convenience, though 
initialization to some other value is straightforward. 
For BURMRAN the array r must be specified. 

Each program uses two random number generators 
employing a machine independent congruence technique (28]. 
One is used to generate uniformly distributed random 
numbers between 0 and 21T for 9 and the other is used to 
generate samples of n It is assumed that 6 15 
Uniformly distributed between + /3 v, thus in the limit 


ав М-оф, E[6 1 = 0 and E[6>] ВА АПОУ САС ОП ОЕ the 





=1 /4= 


program revealed that if the number of paths 15 greater 
than three, phase random multipath properties are 
obtained. Because of the central limit assumptions 
inherent in the phase random model [9] (e.g. N>3), the 
exact nature of the density for Ж is unimportant. 
RANDPHASE was run with Au distributed normally without 
any change in the results. 

The number of independent samples of Des and 8 
required for any given simulation is given by the product 
BEL, N, and I. For most of the simulation runs 
I = 600 samples. M L and N are large, the number 
of independent samples of des and 8 can become quite large. 
One run of BURMRAN simulating the N=15 problem of 
Reference [4] required % 400,000 ПСЕ samples of 
E. and PE The maximum integer that can be accommodated 
on the Interdata 80 is 32,767. I performed a run test 
on the system random number generator (which employs the 
linear congruence method) and for some primes the maximum 
string of unrepeated numbers was 8192, clearly inadequate. 
in order to increase the cycle lengths the machine 
independent congruence method is used. This technique 
artificially increases the maximum allowable integer 
number with factors appropriately segmented so that 


numerical overflows can be computed without machine 
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overflows. The cycle length is greater than 197 and for 
E x 109 samples generated by this method, the 
Chi-square goodness of fit test [29] for uniformity was 
passed at the a = .05 level of significance, though only 
by a small margin. It is true that the leftmost digits 
are the most random and that for large numbers of samples 
some correlation occurs as more of the rightmost dicits 
are repeated. Table IV, reproduced from Reference [28], 
gives the repeat characteristics for the two sequences 

or generators used. Table IV indicates the randomness 

of the generator as a function of cycle length or number 
ОЁ samples. As more and more samples are taken, i.e., 
the cycle length becomes larger, more and more of the 
rightmost digits are repeated which can result in a 
correlation that will introduce some error in the 

overall simulation results. Because the repeat cycle 
length is so large, however, for most Simulation runs 
this error is small. 

When modulation or finite bandwidth effects are 
Simulated, the Interdata 80 random number generator is 
used in addition to the two described above. Because L-1 
for these simulations the large number of Bee 9, 
апа 8 is not required. 


When the programs are executed, histograms are 





Cycle 





TABLE IV 


Sequence One 


735776465527, 


062221556427, 


650107434527, 


617512155527, 


65527, 


465527, 


6465527, 
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Repeat characteristics of the random number 
generater for the two sequences used. 


Sequence Two 


311037552421, 


113326416521, 


655503553421, 


514633562421. 


52421. 


552421. 


1552421, 





с = 


produced by separating the range of samples for each of 
the variables specified by Equations (2.1) and (1.164) 
into 25 equally spaced bins. In RANDPHASE each histogram 
is then plotted along with its respective theoretical 
density and the Chi-square goodness of fit test is 
Bpplied. The area under the theoretical curve is 
computed numerically using the trapezoidal rule, each 

of the twenty five bins being further subdivided into 

16 intervals each. The expected frequency is then equal 
to the number of samples times the area. Bins in the 
tails are grouped such that the minimum expected 
frequency is 7. Thus the number of independent class 
intervals (a function of the variable, number of sources, 
paths, etc.) usually varies between 15 and 25 which is 
close to the criteria specified by Bendat and Piersol [29] 
for % 600 samples. The Chi-square statistic (x?) is 
computed and compared to the pass/fail value Cex S05) 
where n indicates the number of degrees of freedom of 

the statistic and .05 is the level of significance. The 
hypothesis is accepted if y? « e: COS Number or 
degrees of freedom n is obtained by taking the number of 
class intervals and subtracting 1+8 where 8 is the number 
of independent parameters that are varied to fit the pdf 


BO the histogram (e.g. 8=0 for the uniform pdf and 8= 1 
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for the Rayleigh, Gaussian, and Longuet-Higgins pdf's). 
It should be noted that the theoretical densities 


were not least-squares fit to the histograms. Given the 


2 
it 


vis obtained by taking the expected value of 82 From 


input parameters, we know o; from Equation (1.2), and 

the computer generated 8115 (which is always within a 

few percent of the input value of v). The theoretical 

densities are plotted and the samples fall where they may. 
In addition to the output on the graphics display 

a printout for each variable for each run of the simulation 

includes: the variable simulated, the number of sources, 

paths, samples, number of expected and observed frequencies 

for each bin, including an overflow bin, the observed 

and theoretical mean and variance, the range of values 

obtained, and the actual Chi-square statistic and the 


pass/fail value. 


2.2 Simulation Results for a Single Narrowband Source 


Examples of the results of the computer simulation 
RANDPHASE for a single, narrowband, non-modulated source 
are given in Figures 13-36. The figures are plotted by 
variable for N = 4, 5, and 12 paths. For all these runs, 
the input value of v is .007 rad/sec and I = 700 samples. 


Bible у compiles the Chi-square statistics (x^) for each of 
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the figures and the m: .05 in which n is the number of 

degrees of freedom, which is equal to the number of 

class intervals minus one (recall that the densities 

were not fit), and .05 is the level of significance. 

For x? « xi: .05 the theoretical density passes the 

Chi square test. Also shown in Table V is the 

theoretical value of the mean and standard deviation and 

the measured values from the computer generated "data". 
First, note that all the theoretical densities match 

the histograms very well as indicated by results of the 

Chi-square test. Particularly note that the pdf's for 

X and À given by Equations (1.4) and (1.7) respectively 

are indeed supported by computer simulation. Of the 

24 examples shown, only two fail the Chi-square test. 

All the pdf's in Figures 13-36 are the limiting densities 

as N+, The results of the simulation for N=2 and 3 

paths (including further analysis of this case) are 

given in Appendix C. Upon comparison with the results of 

Appendix C, it is quite acceptable to assume that for 

N » 4 paths phase random multipath propagation is 

obtained. It is true that Х for N24 paths did fail 

the Chi-square test but this is only one out of eight 

variables. It should be noted that for N22 and 3 paths 


the limiting densities for seven out of eight of the 
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ariables (the multipath phase 9 is always uniform) 

811ей the Chi-square test (see Reference 232), thus>the 
ustification for choosing N>4 as the required number of 
aths is evident. The failure of ф for the 5 path case 
ndicates perhaps a lack of randomness or a correlation 
or that particular cycle length in the random number 
enerator. Though the results on the average get better 
S N becomes large, as indicated in Table V, (except 
uriously enough БОГ A), after 12 paths the results 
re approximately constant indicating that the 
orrelation in each of the random number generators is 
ШІ -пе !1тплЕјпд factor. This assertion follows from 
Onsideration of N =30 [27] in which by the way x^ for 

ES 4.49, and ae „05 = 26.30. In addition to 
ontaining the other simulations not shown in m figures 
or mentioned above Reference [27] also includes the 

Ore detailed printout mentioned earlier. It is 
nteresting to note that even when phase random conditions 
me met for the single source the pdf for p always appears 
> be among the worst performers while 6 is always among 
ne best. This result was first noted by Hamblen [9] 
men he investigated ocean acoustic data. Because this 
sult is repeated in the simulation, it appears to be 


ENEnction of the statistical nature of the variables, as 








= 130 


well as the non-stationarity of 01 when analyzing real 


data. 


2.3 Simulation Results for Multiple Sources 
Шцшсегз 37-42 are the results of a simulation run, 


ЕП 20 equal intensity, equal v^ receptions, 5 paths, 


and I - 600 samples. The theoretical densities are given 


by Equations (1.19), (1.67), (1.79), (1.80), (1.83), and 


(1.84) respectively. The results of the Chi-square test 


are tabulated in Table VI. Also shown are the means 


апа standard deviations as predicted by the theory and 


observed in the simulation. With one exception, the 
computer simulation agrees well with the theoretical 
predictions. One of the important results of the 
Bnalysis of the equal intensity, equal v^ case 

EE rion 1.2.2.1) is that Рр, (0) is.independent of the 
number of sources. Figures 24, 43, 44, and 40 are 
examples of P. (6) for 1, 3, 4, and 20 equal intensity 
(91 Zora, eaual v* (v = .007 rad/sec) sources. 
Figures 43 and 44 passed the Chi-square test, and the 
histograms clearly support this prediction as well. 


As a final example of multiple source simulations, 


I treated the N=15 group problem of Reference [4]. As 


Was noted in Section 1.2:1.4, and illustrated in Table I, 





Ig. 


B7 
38 
B9 
40 
41 


42 


SS 


TABLE VI 


Results of the chi-square test for RANDPHASE 
simulation of 20 equal intensity (0] = 2.5), 
equal v^ sources with N = 5 for each. The input 
v was .00700. Actual v was .00699 rad/sec. 

The units are the same as in the figures. 


Mean 
Variable х? к 05 
Theory Data Theorv Da 


Std. Dev. 


Ca 


X 9.26 26.30 100.00 99.86 OS 
X 17.01 23.68 0.00 -0.01 0142280192) 
0 11.08 26.30 9.94 9.94 ШОО 5130 
0 19.62 18.31 0.00 0.00 001 001 
А 15.26 25.00 19.89 19.90 0.98 0.93 
Ac 13.53 26.30 0.00 0.00 ШЕ О! 
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the first eight groups are the most energetic. I 

"applied BURMRAN therefore, with 123 sources, the 
appropriate r array, 5 paths, I = 600 samples, and an 
input v of .007 rad/sec. The results are shown in 
Figures 45, 46, 47, and 47a for x, X, A, and D 
respectively. In Table VII I have listed the Chi-square 
results and other pertinent data. The theoretical 
densities are given by the Edgeworth's series approximations 
EU rIons (1.45), (1.108), and (1.53) for x, x, and A 
respectively. For 6 I have plotted the Gaussian using 

the variance obtained from the data (the curve was not fit). 
ЕЕ excellent performance confirms the assertion made 

in Section 1.2 that even for the most complicated 
realizations of Case (c) the pdf for ó will be Gaussian 

ШИ which the variance is in fact = E[ß°]. The results 

for x and A are not very good (once again the rate 
variables out-perform the amplitude variables). The error 
appears to be in underpredicting the mean as indicated 

-in Table VII. However, these results should be inter- 
E їп light of the very large number of samples 
required for 8, and $ (369,000 for each) certainly 
introduced some correlation error due to the repeat 

| characteristics of the random number generator (see 


Table IV). Both e and 3 barely passed the Chi-square 
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TABLE VII 


Results of the chi-square test for BURMRAN 
simulation of 123 different source problen. 
The units are the same as in the figures. 











Mean Std. Dev. 
I Pig. Variable s EE . 05 
Theory Data Theory Data 
45 X 36.46 23.68 | Dec 5. 30 19,63 1.43 
46 X Па 2368 ИШЕ ГТ ОС боот 
47 A 36.34 21.59 ПБК 5 БО 
47а D mo OD 30.14 0.00 aro -- 0-003 
*Егош equ. (1.60a) which is the Gaussian assumption. Note 


‘that these numbers should not be compared with those given 
in Section 1.2.1 because this simulation included only the 
first 8 groups so there will be a slight difference in the 
mean values (see Table I). 
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Fig. 47a Computer generated histogram of p for the 
L :123 different source problem modelling 
distant shipping noise off Bermuda, andthe 
Gaussian using a variance given by E[p°] 
obtained from the computer generated 
samples. 
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test with results of x? = 994.16 and 981.53 respectively, 
and Ше: .05 = 1019.94. Іп апу case, the non-Gaussian 
skew is certainly evident in the histogram in Figure 45, 
once again dissuading the Gaussian assumption in favor 


of the Edgeworth's series. 


2.4 Simulation Results for a Finite Bandwidth 
and/or Modulated Source 


Figures 43-51 are examples of the RANDPHASE 
Simulation of фи when the source is undergoing varying 
degrees of sinusoidal phase modulation. The histograms 
Mee plotted with Equation (1.169) (Note that v and ёс 
are both given in rad/sec). Figures 52-55 are examples of 
varying degrees of uniform frequency modulation, or 
alternatively varying amounts of bandwidth on ды The 
histograms are plotted with Equation (1.176) (Note that 
м апа А аге both given in rad/sec). In Table VIII I 
have Ee led the results of the Chi-square test, and the 
values of fo, A, or B, as appropriate. All these 
Simulations were run with an input value of v = .007 rad/sec, 
5 paths, and I = 700 samples. With the use of Table VIII 
compare Figure 52 (for B=v) and Figure 53 (for B=2v), 
with Figure 28, the non-modulated pdf for d. This shows 
quite clearly again that only for B << 2v can bandwidth 


effects be neglected. As with the other simulation 
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ЕЕЕ ЧАТ 


ЕР <= ос Chi-square test and other data for the 
RANDPHASE simulation of фу for varying degrees of 
modulation and bandwidth with N = 5 paths. 
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48 19.30 2368 ‚0070 = = = ОДО 
49 81.55 23268 .0070 = = .0970 
50 2.13 25268 .0070 = = O 
рі. 22239 DO „0070 = z .0840 
52 25.82 26.530 „0970 01035 0970 = 
Б 3 12.72 23208 „00 70 ‚0070 0140 = 
54 197.62 232688 .0070 5 612 0. 0420 = 
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results, the fit of the theory to the computer data is 
quite good with only one of the eight runs failing the 


-hi-square test. 


An, 
CHAPTER 3 


DATA ANALYSIS 


I have analyzed two sets of data from acoustic 
experiments in the ocean. These data are compared with 
the theory presented in Chapter 1. In Section 3.1, data 
from an experiment performed by R. Porter and R. Spindel, 
near Eleuthera [12] are shown among other successful 
comparisons to support the theoretical pdf's for 
amplitude rate and level rate derived for the first time 
in this thesis. Also, the analytical results of the 
modulation theory derived in Section 1.3 explain heretofore 
unobserved phenomena of the Eleuthera data accurately. 

In Section 3.2, data from the CASE experiment [13] 
conducted in the Pacific are investigated in general, 
with special emphasis in light of the modulation theory 
of Chapter 1. Furthermore, these data provide more 
Ensight into the parameter v^ and reveal, as well, 
shortcomings of our current understanding of this vital 


Ocean acoustic parameter and its driving mechanisms. 


3.1 The Eleuthera Experiment 
Data made available by the Woods Hole Oceanographic 
Institution were acquired during a long-range acoustic 


propagation experiment conducted near Eleuthera [12]. 





E 


The data analyzed consisted of four records (see Table IX) 
in which two cw signals, one at 220 HZ, and one at 406 Hz, 
for each record were transmitted from Eleuthera to 
drifting sonobuoys approximately 300 km northeast towards 
Bermuda. A Doppler position-tracking system [30] was 

used to remove mean multipath phase-rates due to 

sonobuoy motion. These data were also analyzed by Dyer 


and Shepard [15] and Hamblen [9]. 


TABLE IX 


Мес о: the WOH:O.I. Fluctuation Data 


Record Date Time (GMT) Record Length (h) 
447 Ше Sept. 74 1900 - 0310 2216 
448 12 Sept. 74 0402 - 1204 8205 
449 11 Sept. 74 1800 - 0145 STS 


424 ll Sept. 74 0644 - 1527 асло 
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КИ Analysis of the Single Source 


From the original digitized time series of pcosdb and 
popsind, time series of A, X, and $ for the two frequencies 
are produced for each record. The phase data were 
analyzed previously by Dyer and Shepard [15] and 
Hamblen [9]. Details of the phase unwrapping routine 
can be found in Reference [15]. Histograms are then 
generated from these time series. For the histograms of 
À and x, Equations (1.7) and (1.4), respectively, are fit 
Buch that the value of v? and сіу, respectively, minimize 
the Chi-square statistic. The value of v? for each run 
is also obtained by fitting Equation (A8) to the 
histograms of ф. The values of v? obtained from 6 
agree closely, as expected, with those obtained previously 
by Hamblen [9]. However, the value of oe obtained from 
N BpoEEing Equation (1.7) is smaller by a factor of 
almost 2 in every case. This can be en as 
follows. As the analysis of Section 1.3 shows, any 
finite bandwidth effects or carrier instabilities, or in 
this particular case, any errors in the Doppler 
position-tracking system, would manifest themselves in 
Espreading" or increasing the variance of б while leaving 
А unmolested. 


2 


To test this explanation further, v^ is obtained for 





22-5 


2 


each record by taking the ratio of 042 to of, parameters 


| which were found previously [9] by fitting Equation (Al) 
and (A5) to histograms of p and 0, respectively. Also, 
(y? is obtained by taking the value of 21У from the fit 
If Equation (1.4) to the histograms of X and dividing 


by the previously found values of 91 [9]. If the 


reasoning is correct, then these values of v* should agree 


i with the values of v? obtained From A for each record as 


they are all derived from densities which are insensitive 


| to angle modulation or finite bandwidth effects. The 


— и 


results are tabulated in Table X, which shows that 


without exception these values of v^ agree to within 33% 


mor less and, furthermore, they are all as hypothesized 


less than the values of v?^ obtained from ф. 
Except for Record 424, the values for v obtained 
from A, noting the square root range dependence predicted 


by Dyson, Munk, and Zetler [11], are reasonably consistent 


with the values for v obtained previously by them for 


the Eleuthera to Mid-station and the Eleuthera to Bermuda 
transmissions for the 406 Hz source, while the values 

of v obtained from b are not as consistent (see Table XI). 
Why run 424 exhibited the larger values for v in all 

cases (see Table X again) is unknown. 


In order to obtain a measure of the bandwidth of the 





(! 
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TABLE X 


The values of v? obtained from the time series and 
histograms of the amplitude densities are compared 
with the values of v^ obtained in fitting 

Equation (A8) to the histograms of ó. The values 
Of A, the half-bandwidth measure of the 
"modulation" is tabulated for each run in the last 
column. 


eae v? rad?/sec? (x 1077) 
Record BE > 2. 44242 : : 
(Hz) 01% /сү (91/97) A ф A, mHz Io mHz 
220 1.8 T9 240 %7 a3 2.0 
447 
406 3.5 3.6 Wolf 27.9 4.3 2.6 
220 6.7 8.3 6. LORS. 202 киз 
448 
| 406 2.9 254 S 721 4.4 2557 
220 24 2 и 22 9 4 2.6 1726 
449 
406 Sl ZU 3.6 1.6 4.4 2.7 
220 10 1? ЕТ 15 2. TA 
424 


406 11 11 15 24 5.8 326 





ЕСЕ 


TABLE AI 


The average values of v in mHz for the 406Hz 
source for Records 447, 448, and 449 
obtained from 4 and A are compared to the 
values of v obtained by Dyson, Munk, and 
Zetler (DMZ). 











ф A |^. DMZ DMZ 
Mid-Station Bermuda 

Range (km) 300 300 55/0 1250 
Measured Equ. Ray Mix Measured 


v in mHz 4.4 3.0 2.8 (Із) 26 4.0 


242 





[and find o 


c ILS 


Modulation in these data, I fit Equation (1.176) to the 


| histograms of Ф using the values of v?^ obtained from the 
histograms of À. I determine the value of A, the half- 


| bandwidth of the modulation, that minimizes the Chi-square 


Statistic. I likewise fit Equation (1.178) which assumes 


NB Gaussian rather than a uniform modulation function 


G! the standard deviation of the carrier 


fluctuations. Note that the half-bandwidth of the 


modulation as given by l/e of the best fit of 


— 


Equation (1.178) is approximately equal to the value of 


A obtained from fitting Equation MATO The values of 


. v? from the A histograms are used because A, as revealed 


by Equation (1.7), is independent of 93 and is less 
sensitive to its non-stationary behavior. In Table X, 


the least Chi-square value of A and gn are given for each 


G 
run. Because the frequency stability of both sources 

is many orders of magnitude less than v, the mođulation 
as measured by A or сс сап be attributed to the error in 
the Doppler position-tracking system. An error of +.017 
m/sec in measuring the velocity of the receiving 
hydrophone itself results in a half-bandwidth value for 
the modulation of 2.5 mHz at 220 Hz and 4.5 mHz at 406 Hz, 


consistent with the values tabulated in Table X. The 


velocity error is consistent with the experimental setup 
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and the performance characteristics of the Doppler 
position-tracking system [31]. 

In Figures 56-59, the histograms of x, A, and $ апа 
the least squares fit of Equations (1.4), (1.7), and 
(1.176), respectively, are shown for all the records. 
Table XII lists the Chi-square statistics for À and x 
and Table XIII lists the Chi-square statistics for 6 when 
fit by Equation (A8) (which assumes no modulation), 
Equation (1.176), and Equation (1.178). With a few 
exceptions, the fit of the theory to the data is 
excellent. Of the three runs that failed the Chi-square 
test, two were the result of fitting the no-modulation 
density, Equation (A8), to $ which iio tient. GL. the 
foregoing analysis, is suspect from the start. As 
Table XIII shows, while Equation (1.176) performed better 
than Equation (A8) in five out of the eight runs,, 
Equation (1.178) performed better than both Equations (A8) 
and (1.176) for all cases except Run 424 where performance 
was better than Equation (1.176) but slightly worse than 
Equation (A8). In fact, a Gaussian-like error in 
extracting the mean phase rate is characteristic of the 


Doppler position-tracking system [31]. 
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TIE XII 


Results of the Chi-square goodness-of-fit test 
ШЕК апа A to Equations (1.4) and (1.7), 


respectively. 
Freq. * 2 2 Мо? pe 

Record Hz ХУ X. 1. 05 Day Xn:.05 
220 Io NS 3.45 21.03 

447 
406 552 DOMOS 12.816 23.68 
220 1.45 31.41 ПСТ 22m7 

448 
406 абл 19.68 4.51 ZB 
220 12:69 2103 1220 23.68 

449 
406 2.96 21 305 6.44 23.068 
220 25230 зат So Е 29597 

424 


406 6.04 qom 16938 28.87 





= 2206 


TABLE XIII 


Results of the Chi-square goodness-of-fit test 
9 со Equations (A8), (1.176), and (1.178). 


ес. 


Eq. 


Eq. 


(A8): 76 (LAS): 2. 

Record Hz x? y? y? Xp +05 
220 76.24 3.19 DS 25 3L 

447 
406 16.54 4.13 1.19 210205 
220 20.57 ОД 1.0 7.06 23.68 

448 
406 22.95 6.34 BROS 21.03 
220 5.43 35268 2.67 2.03 

449 
406 14.87 p] 02 5267 ое ОЕ 
220 10.45 23.84 12207 26.30 

424 
406 9.78 14.61 12:71 26550 
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peewee Multiple Sources 

Because the 220 Hz and 406 Hz sources were transmitted 
and recorded simultaneously for each record and were each 
quadrature demodulated about their respective center 
frequencies and summed incoherently, this affords an 
opportunity to check the equations derived in Section 1.2.2.2 
for two different sources. For each record, the 220 Hz 
source was designated Source one, and the 406 Hz, 
Source two. The values for 911 and 012 were obtained 

2 


from Hamblen's [9] previous analysis. The values for vi 


2 are obtained from Table X. In 


and v5 in rad*/sec 
Figures 60-63, I have plotted representative results. 
In Figure 60 the histograms for x and p are plotted with 


the following equations: 


E (X) = ————— E SE) аер (= er] (3.1) 
2D 2011 - 2035 11 12 
0 03 о? 
$5 (р) = NU ехр е pee: mac = 2077) 302) 
2D 9117 912 11 12 


Equation (3.1) was first derived by Dyer [4], and 


1 
72 


Equation (3.2) is simply the transformation р= х of 


Equation (3.1). In Table XIV, I have listed the results of 
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Fig. 60 The histograms of X and p for two different sources from Record 
447, plotted with equs (3.1) and (3.2) respectively. 
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Fig. 63 The histogram of A for two different sources 
from Record 447 plotted with equ. (1.123). 
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Be Chi-square test and listed the values of of. and vi 


2 
121 
used. Note that the theoretical pdf's for these plots 


were not fit to the histograms. Also in determining the 


Chi-square statistic, I have approximated the area under 
the theoretical pdf for each class interval by taking the 
value of the pdf at the center of the interval times the 
width of the interval. Figures 61-63 are the histograms 
mor Y, сб, and A, plotted with Equations (1.95), (1.112), 
and ЕЕ 123), respectively. 

In general, the amplitude densities perform rather 
poorly as expected, while the rate densities perform 
exceptionally well supporting the theory in 
Section 1.2.2.2. As noted earlier, the non-stationarity 
in 97 discovered by Hamblen [9] in these data, as well as 
the apparent statistical stability of the rate variables 


accounts for the difference in performance of the 


amplitude and amplitude rate variables. 


2 Crossing Rate Statistics 
Following Dyer and Shepard [15], I have analyzed 


the crossing rate statistics for each of the frequencies 


of the four records of the Eleuthera experiment. Even 


though the "modulation" in these data are very small 


(see Table X), the effect on the crossing rate statistics 








с 


lis Пс. From Equation (1.18), I have [15], 





2 
С(р) = Fo exp [- E (3.23) 
0,727 = 
and [15],. 
G(o.) = 5 (3.4) 


Equation (3.4) assumes no modulation. When the ratio 
Mee squation (3.3) to Equation (3.4) is taken, v cancels, 


me. [15], 


С (р) 


“С 





= (8r)? Р (р) oy (3.5) 


However, if the source is modulated or has a finite 
| > | 
Bandwidth (i.e., B^ 2v), then in fact v does not cancel 


when the ratio of the two is taken: 





or) y V 
О = (8т) < 2р (р) оу (3.6) 
С ф f 
О 
where Va signifies the value of v obtained from the 
amplitude densities, and Уф the value of v obtained from 


fitting Equation (A8) (which assumes no modulation) to the 


histogram of ó of a modulated source. This is tantamount 





po 


EO assuming that the pdf for M is such that Pata тон (1.165) 
yields the Longuet-Higgins density with a larger 


у = 7% > УА” ЈЕ there is no modulation of the source 


mand в << 2v, then v,=v, (as will be shown in Section 3.2). 


| System. The ratio уА/У 


А $ 
For the data of the Eleuthera experiment, however, 


Va Va because Of the error in the Doppler position/tracking 


$ 


Б ССОГ to account for modulation effects in crossing rate 


can be considered a Correction 


| statistics. However, because the error in extracting 


the mean phase rate is more Gaussian than Longuet-Higgins, 
as demonstrated quite convincingly by the performance of 
Equation (1.178) on the histograms of 6, I have also 
computed the correction factor, uSing the analysis of 
E: 93.3, given by Equation (1.209). in Table XV, 

I have compiled the correction factors given by u 
E uPSbation (1.209). I obtain Val by taking the ratio 
Br v from N couv from $ given in Table X.. For 

Fmation (1.209) OG: and v (from A, converted to mHz) are 
likewise obtained from Table X. In Figures 64-67 
reproduced in part from Reference [15] I have plotted 

the new curves applying the Va Yg! and Equation (1.209) 
correction factors. Note that because the correction 
factors are a function of v and parameters of the 


mođulation there is one curve for each frequency. The 
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TABLE XV 


The correction factors to the unmodulated 
value of G(p_)/G(¢_) for the Eleuthera 
experiment, Based 6n the а ОНО 

Bum Table X), to v. (from 6, 

Table X), and a Gaussiañ frequency 
modulation assumption. 


Frequency Val 4 EG: 
220 265 
406 sue 
220 „84 
406 EC 
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220 286 


406 eo 


o 


de 209) 


„41 


„44 


„об 


‚40 


„47 


„42 


58 


· 54 





25 


2:0 


SES 


RECORD 447 


ЕСІ S) o 220 Hz 
х 406 Hz 







406 Hz ( Vs / Và) 


CN 220 Hz (vj) 


o 406 Hz Gaussian 
N ~ ве 
SS SEN 220 Hz Gaussian 
NN RC 
No 
NN IN M 
So 
RS 
АО mn = 
О 0:5 L.O 52 2.0 Eno 3:0 345 40 
Ро 
1/2 
[po » /2] 


Fig.64 Ratio of phase period to amplitude period, versus amplitude, 
for record 447. Equ. (3.5) is the theoretical ratio for a 
non- modulated narrowband source. The correction factors 
уд 1 (Table XW), and Eau. (1.209) (Gaussion modulation) have 
been applied to obtain the corrected curves for 220 Hz and 
406Hz accounting for the modulation. 
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Fig.65 Ratio of phase period to amplitude period, versus amplitude, 
for record 448. Equ. (3.5) is the theoretical ratio for a 
non-modulated narrowband source. The correction factors 
уд 7 (Table ZW), and Equ. (1209) (Gaussian modulation) have 
been applied to obtain the corrected curves for 220Hz and 


406 Hz accounting for the modulation. 
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Fig.66 Ratio of phase period to amplitude period, versus amplitude, 
for record 449. Equ. (3.5) is the theoretical ratio fora 
non-modulated narrowband source. The correction factors 
Vs / Vi, (Table XV ), and Equ. (1.209) ( Gaussian modulation) have 
been applied to obtain the corrected curves for 220 Hz and 
406 Hz accounting for the modulation. 
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Fig.67 Ratio of phase period to amplitude period, versus amplitude, 
for record 424. Equ. (3.5) is the theoretical ratio for a 
non-modulated narrowband source. The correction factors 
AZ (Table XZ), and Eau. (1.209) (Gaussian modulation) have 
been applied to obtain the corrected curves for 220Hz and 406Hz 
accounting for the modulation. 
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bracketed Seno.) > indicates that an average value of 


 6(Ф6) has been used. The range of values results from 


restrictions placed on the maximum value of ¢ aS a 


consequence of phaSe unwrapping routines. As indicated 


in the figures, Po has been normalized such that the 


| Elp] = 2. Details of the above and other aspects of the 


data reduction can be found in Reference [15]. 


The corrected curves are a dramatic improvement 


Ж ИЕР ЕСПЕ Dyer, Shepard theory. Except for run 424 where 


the GausSian performs much better, the Longuet-Higgins 
and Gaussian curves perform about the same with the 


former tending to be slightly high and the latter 


slightly low. The variance in these results only 


| indicates our uncertainty in the exact nature of the 


modulation in these data. The important results are as 
follows: (1) The modulation theory discovered the 
heretofore unnoticed modulation in these data, 


(2) accounted for its effects on the histograms of о, 


апа (3) соггес+1у predicted the crossing rate statistics 


| completely explaining the consistent over prediction of 


the Dyer, Shepard theory. 


3.2 The CASE Experiment 


Data from the CASE experiment [13] were analyzed in 
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order to compare the theory to data in which the 
modulation is known and controlled. Three configurations 
were used: (1) fixed source on a seamount, (2) source 
towed by a surface ship, and (3) source mounted on a 
submersible. The signals were monitored at four widely 
Separated fixed deep water ere at ranges varying 
from 200 to 400 km. Runs from one of the receivers had to 
be subsequently rejected due to an extremely low signal to 
noise ratio. Two carrier frequencies were employed, one 
at 15 and the other at 33 Hz. The signals were frequency 
modulated by a pseudo-random function generator with an 
average period of 107 sec and an overall pattern that 
repeats itself every 640 sec (see Figure 68). The 
bandwidth of the modulation is pre-selected, and for the 
runs analyzed is either 0 (no modulation), 1/8, 1/4, 

Er 1/2 Hz. 

A total of 88 runs of 15 min duration each was 
analyzed. From the filtered digitized complex data, time 
Series and histograms were generated for N and 6. As 
before, Equations (1.7) and (A8) are fit, respectively, 
such that the value of v? obtained minimizes the Chi- 
square statistic. 

Of the 88 runs, 12 were not modulated and the values 


of v obtained from \ and Ф for those runs agree (run by run) 
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to within approximately 10%. In Table XVI, I have listed 
these values as well as the range frequency and receiver 
number for each run. For all of these runs, the source 
was fixed on the seamount. For 9 of these 12 runs, the 
quadrature components exhibit significant nonzero means. 
During the data reduction these means for all the runs 
analyzed regrettably were removed. However, runs 76-77 
had no significant means and the values of v are 
consistent with those from runs 66-68 which are at the 
same frequency. It appears that extracting the means 
even when significant (and this was true only for runs 
66-68, 71-73, and 113-115) is a negligible factor. 
However, runs 113-115 do in fact contain the largest 
three values of v obtained for all the 15 Hz runs and are 
perhaps suspect, though other runs with insignificant 
Means also exhibit values of v approaching and exceeding l. 
For the remaining 76 runs, 24 were modulated with a 
EN bandwidth, 35 with a 1/4 Hz bandwidth, and 17 with 
a 1/2 Hz bandwidth. As expected with modulation, the 
value of v? obtained from fitting Equation (A8) to the 
histograms of ф were larger than those obtained from 
A by large factors depending upon the bandwidth of the 
modulation. Because the modulation pattern (Figure 68) 


is approximately a saw-tooth and the spectrum of the 





Ed 


TABLE XVI 


The values of y obtained from the 12 non- 
modulated runs analysed in the CASE experiment. 
In addition to being non-modulated the source 
was fixed on a seamount. 


EB. E Р > Range (km) RCVR 
# ф A 
66* 33 .094 .104 250 Al 
6т 33 . 353 .354 320 2 
68* 33 .191 .194 450 3 
71* 15 .407 Eu 250 i 
72* 15 A 52 320 2 
73% 15 .741 029 450 3 
76 33 .042 „044 250 1 
77 33 .154 E 320 2 
78 33 „106 .096 450 B 
113* 15 1.349 1.208 250 i 
114* 15 12095 I 320 2 
115* 15 1.279 1.084 450 3 


*These runs showed significant means in the quadrature 
Components. 
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шешасјоп (Figure 69) is approximately uniform, the 
analysis of Section I.B is applied. 

ation (1.176) is fit to the histograms of > 
Msing the value of v^ obtained from fitting A and varying A 
Such that the Chi-square statistic is minimized. The value 
of A is then compared to the actual signal bandwidth set 
by the experimenter. In Figure 70 the error in percent 
1s plotted against the number of runs that exhibit a given 
error. For the 76 runs with modulation, the average 
error in determining the bandwidth of the modulation by 
this method is 8%. Predictive ability appears to be degraded 
Somewhat when the oceanic fluctuations are of the same 
order as the bandwidth of the modulation, i.e., when 
мА the average error was 17%. It should be remembered, 
however, that the modulation is not exactly uniform and, 
therefore, using Equation (1.176) is an approximation 
Eo begin with, and this approximation is worst when v ~ А. 
Finally, I note that the error is consistently negative, 
i.e., the foregoing method underpredicts. The reason 
For this is unknown. 

Performance on the Chi-square test was also quite 
mood. All 88 runs passed at the a= .10 level of 
significance when Equation (1.7) is fitted to the 


histograms of À. In fitting Equation (1.176) to the 








ZH/ 93MO d 


FREQUENCY, HZ 
Fig. 69 The source signal power spectrum of the 


modulation used in the CASE experiment, 


reproduced from Reference (32). 
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histograms of 6 femme OFLuns wlth modulation, only 2 runs 
ERU ror the 12 runs without modulation, only 1 run 
fails when Equation (A8) is fit to the histograms of >. 
The fit of Equation (A8) to the 76 modulated runs is as 
anticipated extremely poor; 29 of the 76 runs failed 

and those that passed did so by a very small margin. 

In Table XVII, I have compiled the values of v 
from O, A, actual A, range, range rate, source 
configuration, frequency, and receiver number for each of 
the 76 modulated runs. Under source N c 
indicates fixed on Een T indicated towed by a 
surface ship, and S indicates mounted on a submersible. 

Figures 71 and 72 are representative of the 
modulated runs when A SEE Figure zu the time 
series and histogram of ћ and the best fit of 
Equation (1.7). Figure 72 includes the time series and 
histogram of Ф and the best fit of Equation (А8) апа 
Equation (1.176). Reference [33] contains. the plots for 
all the runs as well as the FORTRAN listing of the data 
analysis program for the CASE experiment. 

Comparison of v? from the CASE experiment with the 
Dyson, Munk, and Zetler model [11] is unfavorable. The 
values of v^ obtained at 15 Hz are, in fact, greater than 


those obtained at 33 Hz, contrary to the frequency scalinc 
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Fig. 72 The time series and histogram of $ for the same 
run shown in Fig. 7l is plotted along with the best 
fit of equation (A8) which assumes no modulation 
(a) and the best fit of equation (1.176) which 
includes the effect of source induced modulation (b). 
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proposed in their model, and these values are in turn 
several orders of magnitude greater than the values 
predicted by the model (see Table XVIII and Figure 73). 
However, the rms single path phase fluctuation for the 
frequencies and ranges of the CASE experiment as predicted 
by an internal wave model only, appear to be << 21, and if 
so, the model would not be applicable to the CASE data. 

In order to obtain a very rough estimate of applicability, 
Equation (118) (for small 8), Reference [34], has been 
Moeed in Figure 74 Бог <ф > = (27)*, for range in km 

vs. frequency in Hz. That portion of the figure above 
and to the right of the line indicates that phase 
fluctuations due to internal waves are > 27. The dotted 
lines represent the CASE experiment, the circles 

represent the Eleuthera data. Despite the assumptions 
involved, the good results of DMZ with the Eleuthera 

data (Table XI), and the poorer results with CASE, could 
be explained by Figure 74. It is apparent that some other 
mechanism must account for the fully saturated phase 
random process evident in the CASE experiment. Other 
possible mechanisms could include range rate, other ocean 
dynamic phenomena such as Rossby waves or meso-scale 
eddies, tidal currents magnified by bottom interaction 


at the receiver location, or rough scattering effects. 








TABLE XVIII 


Comparison of the average values of v obtained from 
the CASE experiment to the predictions of the DMZ 
internal wave model. Also shown are the variances for 
the values of v obtained from CASE. 






Varliance of 
Case v 
(mHz) 


*These values were computed leaving out runs 71-73 and 113-115 
that exhibited significant non-zero means in the quadrature 
SOmponents. 


E 
These values were computed leaving out runs 66-68 that 
exhibited significant non-zero means in the quadrature 
components. 
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Although no strong correlation was found between v 
and range rate, it is apparent from Figure 75 that the 
Larger values of v were obtained at 0 range rate. I also 
investigated other correlations. I found none between 
м апа range, or time of day (see Figures 76 and 77 
respectively). Because the CASE experiment was run over 
approximately one and one-half months with experimental 
runs occurring on many different days, tidal effects 
cannot be inferred from the times. At the present time, 
the dates of each of the runs are unavailable, making 
tidal checks impossible. 

In Figures 78 and 79 I have plotted the values of 
v vs. receiver for 15 Hz and 33 Hz respectively. Because 
each receiver was monitoring a given experimental run 
simultaneously I have connected the values of v at each 
receiver by lines for each run. The time of day is noted 
to the left of the receiver 1 values. An "N" next to 
the time indicates that there was no range rate, either 
the source was fixed on the seamount or motionless for 
that run. For 15 Hz (Figure 78), receiver 2 consistently 
sees a higher value for v than receiver 1, and receiver 3 
consistently sees a larger value than receiver 1. Between 
receivers 2 and 3 the results are mixed. For 33 Hz 


(Figure 79), receiver 2 consistently sees a v greater than 
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receiver 3. For 8 of the 12 runs, smaller values of v 

are seen at receiver 3 than are seen at receiver 1. 
Between l and 2 the results are mixed. These results 
suggest the possibility of rough scattering or a tidal 
current phenomena which would be receiver dependent. 

Note, however, that the no range rate runs ("N") for 

both frequencies all exhibit the same pattern, with 

low values at receiver l, the highest values at receiver 2, 
and low values again at receiver 3. For all receivers 
range rate groups v at lower values than the no range 

rate runs, the latter having a consistently higher mean. 
Thus, as observed earlier (Figure 75), there does appear 
to be some small correlation with range rate to the extent 
that either there is range А in which case v is 
independent of the amount, or there isn't, in which case 

v exhibits a higher mean and variance. This is more 
graphically illustrated in Figures 80-85 in which I have 
plotted v vs. range rate now separating the values by 


receiver as well as frequency. 
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Fig. 80 A plot of v vs. range rate for the 33 Hz signal 
measured at receiver |. 
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Fig. 82 A plot of v vs. range rate for the 33Hz signal 
measured at receiver 3. 
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DISCUSSION AND CONCLUSIONS 


The major contributions of this thesis are twofold. 
First, on the theoretical front the understanding of 
the statistics of acoustic signals propagated to long 
ranges in the ocean has been significantly advanced by 
the derivation of many new pdf's, particularly for the 
rate variables and joint pdf's of amplitude and amplitude 
rate, as well as the theoretical development of the effects 
of modulation. Second, the application of this new 
theory to ocean acoustic data has revealed new 
understanding of the effects of finite bandwidths and/or 
modulation on the statistics, as well as a clearer picture 
of the limitations of current models relating v? to ЕЕЕ 
or experimental phenomena. 

To be more specific, the completion of the family of 
pdf's for the single source and particularly the derivation 
ӨЕ P: (Å) has permitted direct measure of v° from amplitude 
quantities which are independent of bandwidth and/or 


modulation effects. However, more research is needed to 


2 
1 


amplitude and amplitude rate pdf's. Although су does not 
appear explicitly in P;(1),non-stationarities in 93 will 


possibly have some small effect. The effect, however, will 


understand the effects of romstationarities in of on the 


2 


1 and not to the 


be related to the rate of change of с 





absolute change which would affect those variables in 
which 33 appears explicitly, so that measurement of v? via 
A is the preferred method. 

The Eleuthera data analysis reveals that for 
measuring v* use of the time series and histograms of à 
and Equation (A8) apply only when B<< 2v, A<<v, and the 
Doppler shift due to relative source/receiver motion is 
<< 29. It is true that relative source/receiver motion 
induces Doppler modulation, which is path dependent as 
each path has a different arrival angle at the receiver. 
However, in forming A, or, in fact, any of the amplitude 
or amplitude rate variables, the mean Doppler is removed 
and it is only the path to path differential Doppler that 
remains. For $, however, the mean Doppler as well as the 
differential Doppler contribute, and naive use of 
Equation (A8) on ф will result in error. Even when the 
mean phase-rate is removed, as was done with the Eleuthera 
data, we discover that a new criterion must be met, namely 
that the error in removing the mean phase rate must be 
<< 2). 

Recognition of the A method of measuring v° is alone 
one of the major contributions of this research. However, 
a more thorough analysis is in order to determine the 


effect of the differential Doppler on the pdf's of 
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amplitude, amplitude rate, and phase rate when there is 
relative source/receiver motion. Because arrival angles 
are small for long-range propagation, it is anticipated 
that this effect will be small if not negligible in most 
cases. For the Eleuthera data, the mean drift rate along 
the transmission path was on the average 100 m/hr [12]. 
Assuming a maximum arrival angle of 14° the maximum 
differential Doppler shift is 1x 107 “н2 at 220Hz and 
дык 107 4ң2 at 406Hz, an order of magnitude less than v 
for these data. For the CASE data, the differential 
Doppler was also negligible for all the runs analyzed. 

In addition to range rate, other mechanisms which 
affect the value of v^ must be researched. Analysis of 
the CASE data revealed that some mechanism other than 
internal waves must account for the fully saturated 
phase random nature of the data. My preliminary 
correlations indicate that a receiver related phenomenon 
such as slope-influenced tidal currents, or rough 
Scattering, would be a good place to start. However, other 
mechanisms including Rossby waves, meso-scale eddies, etc. 
cannot be counted out. I must confess to some feeling oz 
uneasiness surrounding the great run to run variance of 
the measured values of v^ in the CASE experiment. Clearly 


a mechanism of some temporally varying nature must be 





EIS 


accountable. This again points toward a tidal phenomena. 
HNEUeuthera data reveal a temporally stable v^, even 
more so than Hamblen [9] discovered, because his 

analysis includes what is now recognized as a Doppler 
error. 

The analysis of effects o£ finite bandwidth and/or 
modulated signals coupled with the new confidence in 
obtaining v^ minus these effects has revealed a 
potentially powerful technique for separating source and 
Ocean effects in the received signal. This technique 
uncovered the effect of the M. ler error in the Eleuthera 
data with great improvement in crossing rate predictions, 
and predicted the bandwidth of the modulation of tne CASE 
data with an average error of 8%. This M е could 
be refined by use of more sensitive statistical tests 
such as the Kolmogorov-Smirnov test. Also the pdf's for 
many different kinds of modulation could be tried, not 
only to determine such parameters of the modulation as 
its bandwidth, but also the nature of the modulation 
itself when it is unknown, by comparing the performance 
of the selected densities on the histograms. 

Aside from the applications of the amplitude 
Statistics to distant shipping noise problems [4,5] the 


analysis of the statistics of multiple source cases is a 





relatively untapped reservoir especially for the rate 
Бр ш=-1с5. In order to improve predictions in the tails 
of the densities whose exact solutions remain unknown, 
puesxchernov bound or "tilted" density [25] approach 
mentioned earlier should be investigated. I believe, 
however, that the analysis presented is complete enough, 
and the problems remaining few enough to warrant 
investigation into the applications of this analysis. The 
coherent source analysis may be a reasonable model for 
determining the statistics for the signal plus ocean 
noise. The multiple source cases analyzed may well be 
applicable to other noise problems. Table III provides 
a comprehensive summary of the state of this analvsis 
at the present time. 

The analysis contained in Appendix C, relying on 
Some results of earlier investigators in fields other 
than acoustics, completes the solution of the statistics 
for the amplitude for N « 3 paths. For N » 4 I have 
shown with the computer simulation that limiting pdf's 
suffice. More research is needed into the exact nature of 
Ра Ва): the pdf for the single path phase. As 
previously noted by Hamblen [9] and as supported by the 
analysis in Appendix C, the rate densities of amplitude and 


multipath phase for N small cannot be found unless Ps, С 
n 
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is known. No investigators to date have related, to mv 
knowledge, any knowledge, experimental or theoretical, 


about РА D aside from the usual Gaussian or uniform 


n 
assumptions as I have done. In Appendix C I have 
solved for 2 when N = 2 paths making these 
assumptions. 
In summary, the theoretical analysis of long range 
acoustic propagation presented in this thesis has been 
wita data. Lrom ocean Acoushe experiments 
supported by both computer simulation and comparison with 
extremely satisfactory results. I hope that the analysis 


presented and the conclusions reached will be of use to 


other researchers in the future. 
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APPENDIX А 


БЕСЕ РК DENSITIES FOR THE 
SINGLE NARROWBAND SOURCE 


Assembled below are the pdf's associated with a single 


narrowband source and their characteristic functions [14. 


Notation: 

Р (х) is the pdf of random variable у 

S e MD) is the characteristic function of y 
u, is the expected value or mean of X 


is the variance of x 


A. Amplitude Densities 


2 
(1) P (o) = Æ expl- 2—] , р> 0 ; Rayleigh (Al) 
3 = 201 


ы ЗЕ = 
Ho = 91 m ‚co = с1(2 772) 


M (w) exp(- —1—) D. 5(-iucj) 


where D(z) is the parabolic cylinder function RA. 











ЕСЕ 








КООР (ух) = E exp[- E m 5 Exponential (A2) 
X 201 201 
- 2 Dr 
ы, = 201, ЈУ 401 
M (w) = mum 
X 1 - i207u 
(3) Р, (4) Lcx ехр [= - E exp (>) ] DES or 
2201 20] 
Log-Rayleigh (A3) 


D 


| = em Zr ee 
нт e[ln 201 үл; SA € 2 


му (ш) = (202) 260г(1 + ієш) 


where y = Euler's constant = .5772... 


T(z) is the gamma function 24. 








nose 


1/27 0<ф<2т 
(4) Pg C0) - УА отат (A4) 
0 Otherwise 
A 

= то, 0с? = т^ 

E 22 5 
М, (ш) = — exp (-iun)sinwr 

ф ШТ 


B. Amplitude Rate Densities 





. 2 
(59) Р.(р) = E exp[- —ÉL——] , -e«p«e ; Gaussian 
e Оло у^ дса 
JE 1 
(A5) 
= 2а ee 
Ha - 0, ©» суу 


E ШЕ 72-2 > 
м, (ш) = exp( > ш Оуу ) 
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КР. (х) = E expí- XL, , -o<xs® ; Laplace (А6) 
Х Ag Digan 
iL їй 
MIO E OO ys 
Hg x 
di 


М. (ш) ee 
X 1 + 4ш%01м2 


(7) Ps (A) 3 > , -œ< co ; Lonquet-Higgins 





(A7) 


Ms (w) 2еу|шік,(2еху|ш|) 


where Ky (2) is the modified Bessel function of order one. 








S 


Це РОГ по modulation only: 


E I» ER E m m: 
2,9 - E 372 ; <ф < о пете сан 5 (А8) 
2[1 + т 
” = t= œo 
A 
Ma (w) = vlol&, (vlo]) 
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APPENDIX B 


THE ERLANG AND RELATED PDF's 


NE tana pai Equation 1.19, Section 1.2.1.1) was 
derived by A. K. Erlang in the early 1900's in connection 
with waiting times in telephone operations. It is the 
probability of the time until the s arrival in a 
Poisson process, or the density of the ne order inter- 
arrival time [35]. Note that the first order 
interarrival time (L-1) is an exponential pdf. L can be 
extended to include noninteger values by merely replacing 
the factorial of the Erlang with the T function and hence 
the name gamma pdf. Although mathematically valid, the 
simple physical interpretation of the Erlang breaks down 
when noninteger values of L are introduced. In the case 
of distant-shipping noise, L is an integer. 

In the language of queing theory, it is interesting 
to note that determination of the pdf of x for Case (c) 
` of Section 1.2.1 in all its generality is the same as the 
zommLıon for the pdf of the moe order interarrival time 
in a renewal process in which the first order interarrival 
times are independent random variables distributed 
pecerding tO Equation (1.19) with different и. апа Li. 

Though originally derived independently, the 


slightly less general chi-square pdf with n degrees of 








SS Е 


freedom can be derived from Equation (1.19) by merely 
letting a = 1/2 and L = n/2 14. 
If we let a = п/р. апа L = m, and make the 


ce in Equation (1.19), we have the 


transformation R 


density of the short time rms pressure: 


mr, 2m-1 - (m/u..) R 


ШШ) EX ERE e (B1) 
Б Г (т) и" 
х 
шесте Of the statistics for y, m = (и 75) 7, which is 


the inverse of the normalized variance of х. 

Equation (Al) was first proposed by M. Nakagami in 

1943 [7] to describe the envelope of long range h.f. 
radio wave propagation undergoing rapid fading and is 
known as the "m" distribution. When m-l, Equation (Bl) 

is the Rayleigh distribution as expected. One can теле | 
that the phase-random model of long range acoustic 
propagation is analogous to the rapid fading of long range 
h.f. propagation. Thus, the connection between the "m" 
and Erlang, or gamma pdf, is more than just a functional 
similarity; they describe in alternative language the same 


process. 


| 
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APPENDIX C 


ОЕ ZWEITES FORN < 3 РАТНО 


The solution for the pdf of the amplitude of a vector 
which is the sum of many vectors added with random phases 
has been a problem of long standing interest. Lord 
Rayleigh was probably the first to investigate this 
problem in two papers, the first published in 1880 [1], 
and the second in 1899 [36]. He was, however, concerned 
with the limit when the number of summed vectors is large, 
and he derived the density of the amplitude for this 
limiting case which bears his name. The first investigators 
to tackle the problem of small N (note that the 
characterization of summed paths as random vectors is 
mathematically identical) was Kluyver (1905) [37], and 
Pearson (1906) [38], tego connectson with mosquito 
migration! 

For arbitrary N and unequal amplitudes, А;, Kluyver 
[37] obtained the general solution in integral form as 
follows: 


со 


N 
Boten | u glup) I J (uAj)du (c1) 


O 


For the problem of sums of independent multipaths, we have 


that A, Sr for all i and thus: 








or 


со 


МР (0) = 0 | u J, fue) (3, (ur) 1" du (С2) 
О 


Pearson [38] obtains solutions of this equation for N=2,3 
Кот (> 4 he obtains solutions in terms of series of 


Bessel functions. 


Solution for N = 2 
Though I derived this solution independently, it was 
Pearson [38] who obtained ШЕШІНЕ гог two 


paths I have for the quadrature components: 


Ба 
Ho 


COSO: + COSO) 
: 2 (сз) 


на 
| 


r(sin®, + 5119.) 


where 9, and 9, are independent random variables distributed 


uniformly between 0 and 27. Forming the sta mean square 


pressure x from Equations (C3), 


y = r^ (cos^0, T 2cos8,cosé. + cos*e, + sin^e, 
+ 2sind, sind, + віп 6.) | (C4) 
Using some trigonometric identities I have 
mir, 2r^cos (6, - 8.) (C5) 


Taking advantage of the symmetry of the cosine the random 


variable defined by 9 -8, behaves as if it were uniformly 


J 


distributed between 0 апа т. Thus, forming the cumulative 


Gesteribution function: 
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Р(Х) = зе Р. (х) (c7) 


ШШШ атп the final solution: - 


1 1 2 
P(X) == —— O Fr (C8) 
where, 
2 2 4 
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н. X Е 


2 1 
Bee ; О <р < 2: (C9) 
2 0 T /Ar2-g2 
Ho = Аг/т, 3 = e - M ae 
and, Me 
D 
P, (A) = {E A 
Бет ларе 


2 2 
uy = 2elnr, 0% = e [nm /6 - 2(%n1)?; (C10) 


Prom Equation (C5) I have 
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It is clear from Equation (Cll) that for small N exact knowledge 
ape ЈА is required before rate densities, or joint 
densities of amplitude and amplitude rate can be found. 

I will now perform the calculations recuired to 
find 2; 0 assuming P6. (0,) is Gaussian and then assuming 


Ра (ду) Poe uniform. "POr the former: 
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C EN ои 5 ee (C12) 
Sn n 5 n 
/ 2ту 
ToS pdaf for t = sin (9, - 8,) is given by 
P,(t) = Zn ; а еш (C13) 
о 
/l-t 
The pdf for z = D = B. is by inspection 
2 
c 
2 
4v 
Р (2) = E e ? |2| < > (С14) 


ШЕ 


Because z and t are indevendent random variables the joint 
density of z and € is the product of Eqns. (C13) and (C14). 
I seek the pdf for the product u - zt, and explorting 
symmetry, 


| Ми 


"RR 5 
Р (ш) = | CE (C15) 
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Differentiating under the integral sign I obtain 





exp (- ) 
Д 1 46202 


= туут ѓо mu 





dt (CIE) 


I now make the change of variable u 7 u* and then integrate 


and make the final change of variables X - 2r^u, to obtain: 





“2 A 
22; 00 = —— exp - 8 m. 


дт туут 
|x] < œ (C17) 


where К (2) is the Modified Bessel function of order 0. 


For 6 uniformly distributed 





1 • 
| | AREE 
Pa (È) = 2/3v n (C18) 


n 0 elsewhere 


Note the factors in Eqn. (C18) insure Е[6 *] = US as required. 


Now the density for z = O, - A is 
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Р„(2) = (C19) 


0 elsewhere 


The pdf for t - sin(9, - 9.) ше given by Fan. (C13) . 
Defining u = zt I integrate over the joint density as 


before. The integration is non trivial and requires 
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Garetul concern for the limits: 
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In each case I perform the integration over z first and then 
differentiate under the remaining integral, and integrate. 


After making the final change of variables x = 2r^u I have: 
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> | Ar“ V3v 
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n Pe 
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| xl. 
= 2 |] - Sage |х| < 44 У3у (C21) 
4r^ 3v 


(5-2 solution тог N = 3 





For this case following Pearson [33], 











0 , elsewhere 


where 3 
2 LOr o 
bom 
(Orr) (3r=p) 


SOS 


(С 22) 


апа (=, a) is the complete elliptic integral of the first 


kind. Note that numerical integration is now required to 


obtain the moments. 


Mransromnine, to x ena A: 
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Т 
——o on БК (=, В), A sr 
2727/2, 1/4 2 

B y 274r 2,1/4 DM 
E , elsewhere 
(22359) 
where 
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(172 + =) 3 (зе ~ „1/2 
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and lastly, 
an 
exp (77 
T 
m5 7372) YK (5, ү), -° < A < 2elnr 
РЕТЕЛ 
3A 
ехр (4є 
РА) = Kk (2, 4 Delnr < К < 2eln3r (C24) 
S 2-т2,3/2 ee 
0 | , elsewhere 
where 
а A 
Е _ or exp (57) 


[exp (5) + 21 [32 - exp (222 ] 


The rate density for x, o, or À assuming a Pa Ge 
n 
are sufficiently more complicated than N - 2 that I have 


not solved for them. 


í.3 Computer Simulation 


Figures Cl-C9 are histograms of the computer 
generated samples of the amplitude variables for N = 2,3, 
and 4 paths. The exact pdf for N = 2, and 3 is shown with 
the limiting pdf. Note that for N = 4 the performance of 


the limiting pdf is good enough (see Table V for Chi-square 
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test results) to warrant use of the lmastungMEkorms for 
IEEE Un Figure C10, I have plotted Equation (C21) with 
the histogram of x for N=2 (recall that the computer 
simulation assumes Pe (6 ) is uniform), as well as the 


n 
n 
ШЕЕ па рағ, Equation (1.4). 
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APPENDIX D 


AMPLITUDE PARAMETER VARIATION 


In this appendix I will briefly introduce the analysis 
involved in amplitude parameter variation. This analysis 
is the first step towards understanding the statistics 
of received signals that are either nonstationary inr 
00 ог when r or 91 is purposely modulated at the 
source. 

For the first case, I consider that r is a random 
variable distributed onl forme between a lower limit of 


EN and an upper limit of a > b. It is evident from 


Equation (2.1c) with L=1 that 


Pay = ко (D1) 
where now p can be considered a random variable, distributed 
Rayleigh according to Equation (Al) but in which 

97 = N/2. Note, however, that there is an implied 
assumption that N > 4. From Equation {Dl1) and the fact 

that r is a measure of the gain of the signal which can be 
controlled at the source, or of the change in pathwise 
signal strength, it is true that r and p are independent 
random variables. This assumes of course that during the 


experiment the number of paths remains constant. Thus, 
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forming the joint density of r and p I obtain the 


emmulative distribution function, 








Pav Pav 
| | 1 Б О o? ? 
р 0 x. ES | | келе етістер 
Pay ^" a-o E 
PAV b 
a 
Pay 
+ — | | о ехр(- 292) 140 (D2) 
(a-b)o] 1 
0 b 
: : г З д 
Performing the integration and taking PI 25 (9717) 
AV AV 
I obtain the final result: 
p О 
Р„ (рд) = —— 1 |кеғ(--АУ-) - кгг(—АУ_) 
AV a- D)o? БО 2 ӘС 2 
1 1 1 
(D3) 
91 fm a+b зт Зар 
y = (atb) —/>,0 = 0 most Sr sn 
Pay 2 2 Pay 2 1/3 2 3(a+b)? 
Note that when a=b=r,y апа а assūme their no 
PAV B 


variation values (see Equation Al, Appendix A), and 
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likewise, 
lim P (ВЕТ) 28 (р) 
a>r Pav AV 0 
b>r 
With the appropriate transformations, I get: 
-% т--- Dm 
(Xay? - XAV Хау 
Р (Xay? = Doms) > ЕЕ ) - E E ) (D4) 
Хау 1 ро; у2 ас. у 2 
1 E 
A 
qe exp (7) Е exp (>) 
P (Av 775 ad 5 \Erf 
Any AV € © a ) bo, /2 
Aay | 
exp (I. 
- Erf | — : (D5) 
a0./2 
1 
Note: ©] = N/2 in Equations (D2-D5). 


In order to compare the effect of certain levels of 
uncertainty in r on the pdf for p, the following 
procedure must be employed. Given that r is uniformly 
distributed between b and a, its average value is a+b/2. 


2 N ,a+b,? 


Pilot Equation (Al) with Cs (бряг and then plot 


Equation (D3). If this is not done, then one may 
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erroneously conclude that the effect of uncertainty in 
r is quite dramatic when in fact it is not! 
I have taken various values of a and b such that 
22 = 1 and plotted Equation (Al) with Equation (D3) in 
Figure Dl. Note that with 50% (3dB) uncertainty іп с, 
the variation from Equation (Al) is not terribly 
significant; however, for uncertainty > 503, the effect 
becomes very noticeable in the pdf. In Figure D2 I nave 
plotted Equation (A3) with Equation (D5). As before 
with 3dB of uncertainty, the effect is small. However, 
for ^ 7dB or greater uncertainty, the effect is large and 
in fact is very close to the limiting form (100$ or 
edB uncertainty). Nakagami [7] obtains a similar result 
in relation to uncertainty of the mean intensity in 
dB of multipath RF propagation. In obtaining the 
limiting pdf (i.e., 100% uncertainty), it is necessary to 
let b=0 and a=2r in Equations (D2-D5) (note, Erf(~) = 1) 
and rewrite the equations analytically first, otherwise 
considerable computer time and/or overflows will result. 
Next, I assume r varies sinusoidally about a mean 
given by a and an amplitude of b. Thus, 


Р (г) i —————— —— , |r-a| « b 
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No voriation r=| Equ. (Al) 
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\ 20% Uncertainty r=!+.2 


\ 5096 Uncertainty r3 12.5 


A ¡00% Uncertainty r=1+| 


Fig. DI Comparison plot of the effect of amplitude 


variation on the pdf for p. Equ (Al) is plotted 
with equ (D3) with various levels of uncertainty. 
N * 5 paths. 
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The cumulative distribution function is given by 
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Integrating, and differentiating with respect to p and 





AV 
simplifying as much as possible, I obtain 
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PAV бау 
Po (Ce 
AV (а- Б) “03 2 (a-b) ?o? 
(САУ s 
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- ті; | exp (- == EN ——— dy 
1 Pay 2 1 vy(b’-a ) + 220,77 - Piy 
(==) 
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Unfortunately, the analytical solution to 
Equation (D7) has not been found. For the special case 


when a=b, however, the analysis can be concluded: 





Р (Pay) L—— (4) | ехр(- ) 
AV 4a%/1 07 16а 202 
О p 
AV AV 
а ааа "av" 1 zac." 25) 
-> 1 -5 18 


where D (z) is the parabolic cylinder function. As a 


check, the integral of Equation (D8) over Pay iom. to © 


more in fact, 1. 


As a final example, I consider the case in which we 


allow 01 to be a random variable and Equation (Al) to be 


the conditional pdf of p given 01. Allowing с to be 
uniformly distributed between b and a, and applying Bayes 


Theorem, I get: 


О О 
med AV _ "AV 2 
Р (олуу) zer |. I exp ( 201/89] (D9) 
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then: 
2 2 
р О О 
EN AV, _ AV 
d Pay’ = А25 m s RA M 
AV 
where E, (2) is the exponential integral. Applying the 
entity [18]: 
со n_n 
(-1) 2 
NR cu EZ ) И 22521 T) 
п=1 
where y is Euler's constant = .5772... I obtain 
Pav E ЈЕ 
а n 
Р (Pap) = == уп = + ) (-1 (—— — 
Pay AV a-b Б Бе СБ ! 
со oe n 
AV I 
- ) Cor) (D11) 
Be ао 
Anota nata mot plot bananas Ehe numbers I 
obtained for relative levels of uncertainty in 0? agree 


1 
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closely with the results of uncertainty in r, and the 
effect of oF randomization (i.e., < 3dB uncertainty) on 
the amplitude densities appears also to be small while 
Мише forms are Obtained for uncertainty > 7ав. 

A more thorough investigation of this phenomena 
could uncover the statistical dependencies on 
identifiable non-stationarities, such as propagation 
loss due to changing ranges or other oceanic or 
experimental factors which result in the temporal 
dependence of r or ВЯ The effects on the statistics 


can then be analyzed using the procedures employed above. 
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